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Abstract 

Euler-Poincare equations are derived for the dynamical folding of charged molec- 
ular strands (such as DNA) modeled as flexible continuous filamentary distri- 
butions of interacting rigid charge conformations. The new feature is that the 
equations of motion for the dynamics of such molecular strands are nonlocal 
C3 when the screened Coulomb interactions, or Lennard- Jones potentials between 

1—1 pairs of charges are included. These nonlocal dynamical equations are derived 

I ^ I in the convective representation of continuum motion by using modified Euler- 

Poincare and Hamilton-Pontryagin variational formulations that illuminate the 
various approaches within the framework of symmetry reduction of Hamilton's 
0> principle for exact geometric rods. In the absence of nonlocal interactions, the 

equations recover the classical Kirchhoff theory of elastic rods in the spatial rep- 
^ resentation. The motion equations in the convective representation are shown 

' to be affine Euler-Poincare equations relative to a certain cocycle. This prop- 

erty relates the geometry of the molecular strands to that of complex fluids. An 
elegant change of variables allows a direct passage from the infinite dimensional 
point of view to the covariant formulation in terms of Lagrange-Poincare equa- 
tions. In another revealing perspective, the convective representation of the 
^ nonlocal equations of molecular strand motion is transformed into quaternionic 

^ form. 
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1 Introduction 
1.1 Physical Setup 

Many long molecules may be understood as strands of individual charged units. Gen- 
erally, the dynamics of such strands of charged units depend both on th.. Generally, 
the dynamics of such strands of charged units depend both on the local elastic defor- 
mations of the strand and the nonlocal (screened electrostatic) interactions of charged 
units across the folds in the molecule. These electrostatic interactions depend on the 
spatial distances and relative orientations between the individual charged units in dif- 
ferent locations along the strand. One important approach to such a complex problem 
is a full dynamical simulation. However, in spite of the importance of this approach 
for determining certain molecular properties, it provides little insight for analytical 
understanding of the dynamics. 

Continuum approaches to the dynamics of molecular strands offer an alternative 
theoretical understanding which is attractive because of the insight achieved in finding 
analytical solutions. Many studies have addressed the elastic dynamics of the charged 
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strands using Kirchhoff's approach [Tj. For a historical review and citations of this 
approach see [2]. Recent advances, especially in the context of helical structures, 
appear in [SI IH El El El IB]- While many important results have been obtained by 
this approach, the generalization of the classical Kirchhoff theory to account for the 
torque caused by the long-range electrostatic interaction of molecules in different 
spatial locations along a flexible strand has not been achieved, although the force due 
to electrostatic interaction has been considered before. See, for example, the article [9] 
which reviews progress in the treatment of charged units distributed along a strand. 
In general, the lack of a consistent continuous model incorporating both torques and 
forces from electrostatic interactions has hampered analytical considerations, see for 
example [5]. 

The present framework does allow treatment of both torques and forces from electro- 
static interactions. We should note that even in the absence of a continuous model for 
nonlocal interactions, it is possible to obtain static solutions using energy minimiza- 
tion techniques. For example, interesting helical static solutions of pressed elastic 
tubes using interactions that prevent self-intersection of the tubes were obtained in 
[TU] . The difficulty in computing the dynamical effects of torque due to long-range 
interactions among the molecular subunits arises because the classical Kirchhoff the- 
ory is formulated in a frame moving with the strand, but deals with a mixture of 
quantities, some measured in the fixed spatial frame and some in the body frame. 
The torque due to long-range interactions then presents a particular difficulty for the 
mixed representations in the Kirchhoff theory, because it is applied at base points of 
a curve that is moving in space. That is, the spatial Euclidean distances and relative 
orientations of the molecules must be reconstructed at each time step during the sin- 
uous motion and twisting of the strand before any self-consistent computation can be 
made of the forces and torques due to long-range electrostatic interactions. 

In fact, even when electrostatic forces are not involved, the motion of realistic curves 
in space is inherently nonlocal because of the requirement that the curve not cross 
itself during the dynamics. In the purely elastic Kirchhoff approach, such nonlocal 
considerations are neglected. Physically, however, self-intersections are prevented by 
the existence of a short-range potential (e.g., Lennard- Jones potential) that produces 
highly repulsive forces when two points along the curve approach each other. 

This paper casts the problem of strand dynamics for an arbitrary intermolecular 
potential into the convective representation introduced in [TT] and applied in 
exact geometric rod theory in [12]. Its methods are also applicable to the consideration 
of Lennard- Jones potentials and the constrained motion of non-self-interacting curves. 

If the curve were constrained to be rigidly fixed in space, and the attached molecules 
on this curve were allowed to simply rotate freely at each position, the theory of 
motion based on nonlocal interaction between different molecules would be more 
straightforward. Of particular interest here is the work [13] where a single charge is 
attached at each point along the fixed filament by a rigid rod of constant length that 
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was allowed to rotate in a transverse plane. These charges were allowed to interact 
locally with other nearby charges that were similarly attached to planar pendula of 
constant length mounted transversely to the fixed filament. 

This constrained motion can be generalized to allow fiexible motion of the strand 
(time-dependent bend, twist and writhe) while also including the degrees of freedom 
of molecular orientation excited during the process of, say, DNA folding. Accord- 
ing to this class of models, a DNA molecule is represented as a fiexible filament or 
strand, along which are attached various different types of rigid conformations of 
sub-molecules that may swivel relative to each other in three dimensions under their 
mutual interactions. The flexibility of the fllament arises physically because the elec- 
trostatic interaction between any pair of these rigid conformations either along the 
filament or across from one loop to another of its folds is much weaker than the in- 
ternal interactions that maintain the shape of an individual charged conformation. 

This paper considers rigid charge conformations (RCCs) mounted along a fiexible fila- 
ment. The RCCs are more complex than the planar pendula considered in [13]. They 
are allowed to interact with each other via a nonlocal {e.g. screened electrostatic, or 
Lennard- Jones) potential. Our investigation is based on the geometrically exact 
rod theory of Simo et al. [12], which is expressed in the convective representa- 
tion of continuum mechanics. The rotations of rigid charge conformations along the 
fiexible filament are illustrated in Figure [T] 




k=3 



Figure 1: Rigid conformations of charges are distributed along a curve. Note that 
this is a spatial representation of the orientation. 

These rigid conformations of multiple charges are allowed to interact via an effective 
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many-body potential representing their screened electrostatic interactions. The non- 
local interactions among these RCCs depend on their spatial separations and relative 
orientations, which are both allowed to evolve dynamically. Thus, the inertial mo- 
tion of a pair of RCCs mounted at any two spatial points r{s,t) and r{s',t) along 
the filament is governed by an effective potential interaction energy that depends on 
their separation and relative orientation. The filament is taken to be one-dimensional, 
although the orientations of the rigid charged conformations mounted along it may 
be three-dimensional. A practical example to which our filament approach applies is 
the vinylidene fluoride (VDF) oligomer [H], which may be approximated by a strand 
carrying a dipole moment whose orientation is perpendicular to the axis of the strand. 
The VDF oligomer strand is straight for small lengths, but it forms complex shapes 
due to electrostatic interactions for longer lengths. 

The theory presented here generalizes directly to the case when the dimension of 
the underlying manifold (filament) is greater than unity, and so is applicable to such 
problems as the motion of charged sheets, or charged elastically deformable media. 
While we present part of the relevant geometry here, we leave its applications in 
higher dimensions for a later publication. 



Plan. The remainder of the paper is organized as follows. Subsection L2 outlines 
the content of the paper in mathematical terms by giving an overview of the various 
spatial representations of filament dynamics discussed here from the canonical and 
covariant point of views. Subsection 1.3 connects our results to the earlier literature. 
Paragraph 1.3.1 relates the theory presented here to the classical elastic rod approach 
pioneered by Kirchhoff. The need to keep track of spatial separations in long-range 
electrostatic interactions requires that we write the dynamics in either the spatial or 
convective representations, as opposed to the Kirchhoff mixed representation. Para- 
graph 1. 3. 2| considers the simplified case when the orientations of the RCCs along the 
curve may depend on time, but the position of any point s along the curve is fixed, 
thereby connecting to earlier work in [13] . Section [2] incorporates the flexible motion 
of the filament into the dynamics by using the geometrically exact rod theory of Simo 
et al. [T2]. The equations of motion are derived in convective form by using a mod- 
ified Hamilton-Pontryagin and modified Euler-Poincare approach to allow 
for nonlocal interactions, in Subsections 3]T] and 3^, respectively. These equations 
are formulated as conservation laws along the filament in Section |4] and their affine 
Lie-Poisson Hamiltonian structure is elucidated in Section [5j Section [6] explains the 
background for the affine Euler-Poincare and affine Hamilton-Pontryagin 
approaches and applies this framework to the dynamics of charged strands. Sec- 
tion [7] introduces a remarkable change of coordinates that decouples the equations 
into their horizontal and vertical parts. Section |8] explains the geometric struc- 
ture of this coordinate change and leads to the covariant Lagrange- Poincare 
formulation. Subsection 8.10 and Section [9] discuss generalizations of the molecular 
strand to higher dimensions. In Section |9} the equations of motion are obtained by 
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an alternative covariant Lagrange-Poincare approach. Section 10 provides a useful 
representation of the convective frame dynamics of the flexible strand using quater- 
nions. Section 11 briefly summarizes our conclusions and sets out possible directions 
for further studies. 



1.2 Mathematical Setup 

1.2.1 Description of the variables involved 

In the Lagrangian representation, the motion is described by the variables A(s,t) G 
5*0(3) and r{s,t) e M.^. The vector r{s,t) is the spatial position of the filament and 
the variable A{s,t) denotes the rotation of the RCC at the point s along the filament 
at time t. By taking the time and space derivatives, we find the material velocity 
{A{s,t),r{s,t)) and the angular and linear deformation gradients {A'{s,t),r'{s,t)), 
respectively. Given A and r, we define notation for the following reduced variables 



Remark 1.1 (Notation) Quantities defined using derivatives in s are denoted using 
capital Greek letters, whereas lower-case Greek letters (except for p) denote quantities 
whose definitions involve derivatives with respect to time. Bold letters, for example 
r, denote vectors in whereas i7 is a 3 x 3 skew-symmetric matrix in the Lie algebra 
so(3). 

Definition 1.2 The hat map : (M^, x) — > (so(3), [,]) is the Lie algebra isomor- 
phism given by uv = u x v for all v G M^. 

Thus, in an orthonormal basis of and w G M^, the 3x3 antisymmetric matrix 
u := u E so{3) has entries 

Ujk = {u)jk = -^jkiu^ ■ (1-2) 

Here Cjm with j, k, I G {1, 2, 3} is the totally antisymmetric tensor density with €123 = 
-|-1 that defines the cross product of vectors in M^. In what follows, we shall employ 
this notation by writing i7 := f2 and u := Q. 





= A- 


'A' 


G so(3) 




= A- 




G so(3) 


r 


= A- 


-v 


G M^ 


7 


= A- 




G M^ 


p 


= A- 


-V 


G M^ 



The physical interpretation of the variables ( 1.1 ) is as follows: The variable p(s, t) rep- 



resents the position of the filament in space as viewed by an observer who rotates with 
the RCC at {s,t). The variables (f2(s, t), r(s, t)) describe the deformation gradients 
as viewed by an observer who rotates with the RCC. The variables (a;(s, t), 7(5, t)) 
describe the body angular velocity and the linear velocity as viewed by an observer 
who rotates with the RCC. 
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1.2.2 The canonical point of view 

The canonical viewpoint of continuum dynamics derives the equations of motion by 
applying a process of reduction by symmetry to a cotangent bundle T*Q endowed with 
a canonical symplectic form. This approach has been extensively studied for fluids, 
see for example [15] for the Hamiltonian description and [16] for the Lagrangian side. 
In hydrodynamics, Q is the product of a Lie group G and a representation space 
V on which the group acts linearly as G x V —>■ V. The dual space, V*, is the 
space of linearly advected quantities such as the mass density or the magnetic fleld. 
The associated process of reduction by symmetry under the action of G is called 
Lie-Poisson reduction for semidirect products. For such systems (in the left 
version), we have the relations 

^ = ^"1^' (1.3) 

where g{t) EG is the Lagrangian motion, ^(t) is the convection velocity, and a(t) G V* 
is the evolution of the advected quantity for a given initial condition oq. Note that 
a{t) is also a convected quantity. For the molecular strand we have g = (A, r) and 
a = (f2,r,p). However, the relations (1.1) cannot be recovered from ( |1.3 ) because 
the variables (f2,r,p) are not linearly advected. Thus, a generalization of (1.3) is 
needed, in which G acts on a by an affine action. Such a generalization is given 
by the process of affine Euler-Poincare reduction developed in the context of 
complex fluids in [17J. This theory, which we recall in Section[6| produces the relations 

^ = , 1, (1.4) 

a = g ao + c{g ), 

where c is a group one-cocyclej^ If we take = 0, then the advected quantity evolves 
as 



Remarkably, the evolution of (f2, F, p) in ( |1.1 ) is precisely of this form for a well chosen 



cocycle. The variables (A(s, t), t(s, t)) are interpreted as time- dependent curves in the 
infinite dimensional Lie group 

G = J^{I,SE{3)) 
of all 5'i?(3)-valued smooth functions on /. 
Remark 1.3 The variables 

(u;,7) = (A,r)-i(A,f), 



^ That is, c satisfies the property c{fg) — c{f) + fc{g), where / acts on c{g) by a left represen- 
tation, as discussed in Section [6] 



Ellis et al. 



Dynamics of charged molecular strands 



9 



their associated momenta 

and the afiine advected variables {fl,r,p) are all convective quantities, see [18]. 
In this context, convective quantities are also called body quantities, since they are 
defined in a frame following the motion of the molecular strand. 
In contrast, the variables 



.<«p<-'):=Ad;,,,,_.(|-,^) (1.5) 



^{A,r) 

are spatial quantities, i.e., they are defined at fixed points in Euclidean space. 



1.2.3 The covariant point of view 



The covariant point of view interprets the Lagrangian variables {A{s,t),r{s,t)) as a 
section a of a trivial fiber bundle over spacetime. The bundle is given by 



X 



7rxp:P = Xx SE{3) -> X, X = [0, L] x M 9 {s, t) 
and the section a is naturally defined by 

(t{x) = (x, A(x), r(x)). 
Carrying out a covariant reduction for the first jet extension j^a recovers the relations 



(1.1) in a natural way. 



From the covariant point of view, the main results and relationships among the various 
sections of the paper may be understood by tracing through the following diagram: 



A{s,t) e S0{3) 



(A,r) is,t)eSEi3) 



SO {3) 



P 



T^XP 



X 



id 



p{s,t) e 



X 



This diagram lays out the sequences of coordinates and manifolds used here to rep- 
resent the molecular strand dynamics. The definitions of the various projections are 
given in Section [7] and Section [8j For the moment we shall only introduce the variables 
and use the diagram to describe the relationships between the two main convective 
representations of filament dynamics used here. 



Recall from ^ 1.2.1 that A(s,t) G S0{3) denotes the rotation of the RCC at a point s 
along the filament at time t. The vector r{s,t) is the spatial position of the filament 
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at and p{s,t) is the position of the filament in space as viewed by an observer 
who rotates with the RCC at The left hand vertical sequence describes the 

fibers over spacetime of the right hand vertical sequence. Therefore we regard the 
variables introduced above as local descriptions of sections of the spaces in the right 
hand vertical sequence. Since the bundles involved are trivial, these local sections are 
sufficient to describe the global sections. For example, p{s,t) G M'^ corresponds to 
the section 

{s,t,p{s,t)) G S = X X M^ 

This characterizes the space S. These sections naturally lead to the use of linear 
maps in place of tangent vectors. Thus, in place of tangent bundles, we introduce 
jet bundles, as described in Section [8] These jet bundles have local representations 
in terms of fiber linear maps from TX into vector spaces covering the section p{s, t) 
in the vector bundles TP/ S0{3), TS, and adP. The space of all such linear maps 
is denoted, for example, L{TX,TTi). With this in mind we can express the reduced 
configuration spaces in terms of the following diagram: 



(p, nds + wdt)5o(3) e L{TX, (idSE{3)) L{TX, adP) 

1 1 

(p, Yds + jdt, nds + ujdt) e L{TX, TSE{3) /SO{3)) 

T7rsB(3)/SO(3) Ttisp/SO(,3) 

(p, p,ds + ptdt) e L{TX, TR^) y jiE 



J^P/S0{3) ^"^^/^o(^), L{TX,TX) 



id 

L{TX, TX) 



where Tnxp/SO^S) and Th-ep/SO^S) denote the tangent maps Tttxp and Ttt^p 
acting by composition on linear maps with values in the fibers of TP/ S0{3). 
This diagram arises in Section [7], where we find that the equations of motion can be 
drastically simplified by a change of variables that passes from the upper to the lower 
horizontal sequence in the diagram. The geometrical significance of this change of 
variables is made more precise in Section [sj In particular, the variables (p, Pt, Ps) 
introduced in Section [t] are coordinates in the J^S component, where E := P/S0{3). 
The variables (r2,a;) describe the L{TX, a.dP) component, where adP denotes the 
adjoint bundle associated with P. The coordinates on L(TX, TP/ S0{3)) recover the 



definitions (1.1) by covariant reduction. Within this framework, we develop a new 



formulation of molecular strand dynamics. The modified Euler-Poincare argument for 
a Lagrangian l{u}, 7, fl, F, p) and related arguments that use the same set of variables 
(p,r,7,ri,cc;) are formulated as sections of TP/S0{3) = TX x TSE{3)/SO{3), 
where X := [0, L] x M 9 (s, t) is the spacetime. 



Remark 1.4 (Defining the convective representation) Formulas (1.1) for the 



variables in the upper horizontal sequence in the diagram above define the convective 
representation of the exact geometric rod theory [12] . We shall see in a moment how 
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these variables in the convective representation are related to the Kirchhoff variables 
for strand dynamics. 

In Section [8] we introduce the structure described by the left hand vertical sequence. 
This allows us to construct the lower horizontal sequence and the right hand vertical 
sequence. When we formulate the problem in the new geometry in Section [8] using 
the covariant Lagrange-P oincare approach, we recover the change of variables 
introduced in Section [7j In the new geometry, the following coordinates are used, cf. 
the definitions in equation (1.1): 





= A- 


-'A' 


Gso(3), 




= A- 


-'A 


G5o(3), 


Ps 


= A- 


-v 




Pt 


= A- 




- A-^Ap G , 


p 


= A' 


'V 


G M^ 



;i.6) 



This covariant Lagrange-Poincare approach is generalized in Subsection 8.10 and 
Section [9] to consider higher dimensional problems such as the molecular sheet, as 
well as problems such as the spin chain that have different microstructures. 



1.3 Connection to previous studies 

1.3.1 Purely elastic motion and Kirchhoff equations for elastic rod 

The results of this paper may be compared to the classical Kirchhoff theory of the 
purely elastic rod, particularly in terms of the available conservation laws. This 
comparison was presented for the purely elastic case, i.e., the Lagrangian / is an 
explicit (local) function of the variables / = /(w, 7, f2, F, p) in |T2|. This work is 
extended here to the case of nonlocal interaction. 

Of particular interest to us are the balance laws for angular and linear momenta. 
For this comparison, we shall use the notation of [H]. For simplicity, we assume that 
the position r{s) along the filament is given by the arc length s. This assumption 
conveniently avoids extra factors of |r(s)| in the expressions. We shall also mention 
here that in order to connect to the Kirchhoff theory, we need to make an explicit 
choice of A(s) G 5*0(3) as a transformation matrix from the fixed orthonormal basis 
{El, E2, E3} of M.^ to the orthonormal basis of directors {di(s), d2(s), d3(s)} describ- 
ing the orientation of the filament (see Figure that IS, 

d,{s) = Al{s)Ek, t = 1,2,3. (1.7) 

There is some ambiguity in the choice of the basis {di(s), d2(s), d3(s)} at every given 
point. The most popular selection of the basis is governed by the so-called natural 
frame. We shall not go into the details of this right now and refer the reader to [9] 
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for a more complete discussion. In principle, we need not have taken this particular 
choice of A, since for rigid charge conformations (RCC), the relative configuration of 
charges is not changed under the dynamics, and the configuration of an RCC state 
at each point s is completely described by a pair (A,r) G SE{3). Taking A to be 
a different presentation of RCC would lead to a transformation A(s,t) —>■ AA{s,t) 
where A G 5*0(3) is a fixed matrix. While our description is equivalent in this case, 
the explicit relation to Kirchhoff formulas is cumbersome. 

We shall note that if the charge conformations were allowed to deform, then A would 
no longer be an element of 5*0(3). Instead, the charge conformation would be de- 
scribed by a general matrix A and a vector r G M^. No explicit relation to Kirchhoff 's 
formulas is possible in this case. 



As mentioned in Section lA^. Kirchhoff 's approach does not allow for a simple com- 
putation of Euclidian distances between the charges unless the spatial length-scale 
of the rigid charge conformations (RCCs) holding the charges at given point t]j^{s) 
is negligible. It is interesting that in the more complex case considered here, the 
equations become formally equivalent to Kirchhoff's equations, provided the effects 
of non-locality are computed appropriately. In particular, one requires an appropriate 
mapping from the convective representation to the Kirchhoff representation, as well 
as some identities connecting nonlocal contributions to the total derivatives of the 
Lagrangian. This mapping is amplified in more detail in Subsection |4] below. 
The linear momentum density p is defined as p(s) = pd{s)r{s), where pd{s) is the 
local mass density of the rod. In that case, the kinetic energy due to linear motion 
Klin is given by 

Kun = \j pMWKsW^s = IJ pd{s)\\A-'r{s)rds = \ j p,(.)||7(.)f d. ■ 

Consequently, the variable p and the linear momentum dKun/ 6^ are related by 

P = = Ap,7 = • (1.8) 

After these preliminaries, we are ready for a detailed comparison with Kirchhoff's 
theory. A point on a rod in Kirchhoff's theory is parameterized by the distance 
r{s, t) measured from a fixed point in space. The z-th component of the local angular 
momentum in the body frame {di(s), d2(s), d3(s)} is defined by 7r*(s) := I*(s)<-J-'(s), 
where u:^{s) is the j-th component of body angular velocity given by t3(s) := uj{s) = 
A(s)~^A(s), and Ij(s) is the local value of the inertia tensor. Note that the inertia 
tensor I(s) expressed in body coordinates is time-independent. Thus the local kinetic 
energy due to rotation is given by 



and hence 

TT = la; = — r — . 
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To write the conservation laws, we need to express the angular momentum in the 
fixed spatial frame {Ei, E2, E3}. To distinguish it from tt which was expressed in the 
body frame {di(s), d2(s), d3(s)}, we shall denote the same vector in the fixed spatial 
frame {Ei,E2,E3} by tz^^I The same convention will be used for all other vectors. 
Thus, ( 1.7[ ) yields 

7r(s) = n\s)d,is) = ms)u;\s)d,is) = I}(s)a;^'(s)Af (s)Efc = n^^^'\s)Ek, 

so the k-th component of the spatial angular momentum is expressed in terms of the 
local body quantities and u)''{s) as 



^(E),fc _ ^feji i _ j^j^i 



A 



6K, 



rot 



;i.9) 



Thus, the vector 7r^^)(s) of body angular momentum expressed in the spatial frame 
is connected to the local body quantities as 



7r(^) = AIu; = A 



6K 



rot 



5uj 



;i.io) 



Remark 1.5 The vector tt*-^) and all other vectors with the subscript (E) do not 
have the physical meaning of the angular momentum in the fixed frame. The true 
angular and linear momenta in the spatial frame will be denoted (see immediately 
below) with the superscript {S). The quantities with the superscript (E) are just the 
transformations of vectors with respect to rotation of the base frame. Hopefully, no 
confusion should arise over this distinction. 

In general, it is assumed for physical reasons, that the Lagrangian in Kirchhoff's 
formulation has the form 



/(a;, 7, 17, r) = Kun{l) + Krot{u:) - E{n, T) , (1.11) 

where E{n,r) is a certain explicit function of f2 and F (not necessarily quadratic). 
In this case, the body forces n = 6l/6r and torques m = 61 /6ft are connected to the 
transformed quantities n^^^m'^^) in Kirchhoff's theory as 



Next, we use formula (1.5) to transfer to spatial frame. Identifying elements of se(3)* 
with pairs of vectors {/j,, rj) G M^, produces a useful formula for the coadjoint action 

Ad^^_,)-i (^, r]) = {An + rx At], At]) . (1.13) 
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Thus, the spatial momenta - denoted by a superscript (S) - become 

(7r(^),p(-)): = Ad^,..)-.(^^ 

= (7r(E) + rxp(^),p(E)) 



5lj' 5j J V Siv 



^61 , SI 61 

A— + r X A—, A— 
07 07 



;i.i4) 



upon using (1.8) and (1.10). Analogously, using ( |1.12 ), the spatial torques m'^"^) and 
forces n^"^^ are 



= (m(E) + rxn(E),n(E)). 



61 , SI ^ SI 

A— + r X A — , A — 
Sn ST' ST 



The conservation laws in Kirchhoff theory may now be written as 

|(7r(^),p(^)) + ^(m(^),n(^)) = (T,f), 



;i.i5) 



;i.i6) 



where T and f are external torques and forces, respectively. Equations (1.16) give. 



componentwise, the following linear and angular momentum conservation laws (cf. 
equations (2.5.5) and (2.5.7) of [9]) 



|p<- + |K.-f)^o, 

I- (7r(E) + r X p(E)) + I- (m(E) + r x n^^^ - L) = , 
at OS 

where F and L are defined as the indefinite integrals, 

/S PS 
f(g)dg and L = / [r(g) x f (g) + T(g)] dg . 



;i.l7) 
:i.l8) 



Opening the brackets in (1.17) and (1.18) gives the balances of linear and angular 



momenta in Kirchhoff's approach (cf. eqs. (2.3.5) and (2.3.6) of [9]) 

5p(E) 



dt ^ ds 

a7r(E) 9m(E) Or 

+ — + 



n(^) = T. 



1.19) 

;i.2o) 



dt ds ds 

To see how these Kirchhoff balance laws look in our representation, one may substitute 



relations (1.14) and (1.15) into (1.16) to obtain: 
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Assume now that the Lagrangian / depends exphcitly on the additional variable p = 
A^^r. This corresponds to potential forces exerting forces and torques. As shown in 
Section |4], in our representation the external torques T and forces f are given by 



(T,f) 



61 



61 



— X p, — 
6p 6p 



1.221 



By using formula (1.13), relationship (1.22) simplifies to 
61 



Ad 



(A,r)- 



6p ^' 6p 



A 



6l_ 
6p 



X P 



A- 



61 
'6p 



X (Ap) 



+ r X A— ,A— 
6p 6p 

r X A — , A — 
6p 6p 



0,A 



61 

Tp 



;i.23) 



upon remembering that Ap = r. 

Remark 1.6 (Potential external forces produce no net torque) 



Equation (1.23) implies that potential external forces produce no net torque on the 



strand. Hence, the nonzero torques T in (1.21) must arise from non-potential forces. 



The conservation law (1.21) is formally equivalent to the classical expressions in 



(1.17) and (1.18), even if nonlocal interaction is present. This equivalence shows how 



the classical results (1.17) and (1.18) generalize for the case of nonlocal orientation- 



dependent interactions. Clearly, the conservation laws are simpler in the Kirchhoff 
representation. However, if nonlocal interactions are present (called self-interaction 
forces in |9]), the computation of the required time-dependent Euclidian distances in 
the interaction energy becomes problematic in the classical Kirchhoff approach. As 
we shall see below in Section |4| these conservation laws may be obtained, even when 
nonlocal interactions are present. Also in Section |4| we show that the nonlocal forces 
are included in the conservation law (1.21) and are expressed in the same form as a 
purely elastic conservation law. 



The balance laws (1.17) and (1.18) are much simpler in appearance than the ex- 



pressions in (1.21), as they do not involve computing (A, r) at each instant in time 
and point in space. Thus for elastic rods, in the absence of nonlocal interactions, 
the Kirchhoff mixed (convective-spatial) representation appears simpler than either 
the convective or spatial representations. However, the presence of nonlocal terms 
summons the more general convective approach. 

Remark 1.7 (Reduction of static equations of motion to the heavy top) 

A famous analogy exists between the stationary shapes of an elastic filament and the 
equations of motion of a heavy top jl9l EH] . In our formulation, this analogy appears 
naturally. This shows the advantage of using the geometric approach, even in the 
study of classical problems of filament dynamics. This paper focuses, however, on 
the derivations and geometric structures underlying the dynamical equations, rather 
than on the solutions of the equations. 
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1.3.2 Reductions for a fixed filament 

We may briefly apply the ideas of the present paper to the particular case of a fixed 
filament, in order to compare the motion equations with those arising in [T3] . 
The analysis of filament dynamics induced by nonlocal interactions simplifies in the 
case when the position of the filament is fixed as r{s) and does not depend on time. 
For simplicity, we shall assume that the filament is straight and s is the arc length, 
so that r{s) = (s, 0, 0)^. The following reduced Lagrangian is invariant under the left 
action of the Lie group 5*0(3): 

1= ^-J oj{s)-I{sMs)ds -^-J f{n{s))ds -^-JJ U{p{s),as,s'))dsds' . 

V ' ^ V ' " v ' 

Kinetic energy Elastic energy Potential energy 

(1.24) 

A nonlocal interaction term appears in the potential energy of relative orientation in 
this Lagrangian. This term involves a variable 

as,s')=A-\s)Ais')eS0i3), 

which defines the relative orientation of rigid charge conformations at two different 
points in space. The variable ^(s, s') G S0{3) is invariant with respect to simultaneous 
rotations of the coordinate frames for s and s', but it is not an element of a Lie 
algebra. In particular, ^(s, s') is not a vector. The presence of nonlocal interactions 
introduces dependence on relative orientation and thereby produces new types of 
nonlocal terms in the corresponding Euler-Poincare dynamics obtained in applying 
reduction by 5*0(3) symmetry to Hamilton's principle. 

Remark 1.8 (Aim of the paper) 

In this paper, the infiuence of non-locality due to electrostatic forces on rod me- 
chanics is studied by using various approaches, including the Euler-Poincare varia- 
tional method. This variational approach leads to an equivalent Lie-Poisson Hamil- 



tonian formulation of the new equations appearing below in (3.3), (3.4). Applying 
the Ad^^^^_i transformation from convective to spatial variables in these equations 
streamlines the form and exposes the meaning of the interplay among their various 
local and nonlocal terms, relative to the Kirchhoff theory. 



Euler-Poincare dynamics Euler-Poincare dynamics for the angular dynamics on 
a fixed filament follows from stationarity of the left invariant total action 




dt. 



Note that this case does not require computation of the evolution equation for 7 
since the filament is assumed to be fixed in space, i.e., 7 = A"^^ = 0. The variational 
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derivative SS for such a Lagrangian is computed as 
SS= [ 5liu;,A,n)dt = 



6I_ 



6fl) dt. 



1.25) 



for the notation S = A ^6 A. As we will see in Section 3.2.1, these variations are 
related by 

6lj = T, + [lj,T] = T, + ad^J2 , 

= S' + = S' + adoS, 

6p = — S X p . 



Substituting these formulas into (1.25) then integrating by parts in the time t and 
one-dimensional coordinate s along the fiber yields 



6S 



61 dt 



d_5l 
dt 5uj 



51 _ d 51 

' 5lo ds 5Vl 



51 



+ a-d*,- — + ad^ — 



dU_ 



s^s')e{s,s')+i{s,s') 



dU , 



X p) ds' , S ) dt . 



51 
5p 

Thus, Hamilton's principle 5S = implies the Euler-Poincare equations, 

51 



1.26) 



d_5I^ A* — 



5p 



X p ds' 



+ / (^-^(^,^'r(^,^')+e(^,^')(^(^,^0) jds'. (1.27) 

Note that these Euler-Poincare equations are nonlocal. That is, they are integral- 
partial-differential equations. 

Reformulating (1.27) in terms of vectors yields the following generalization of equa- 
tions considered by ^3], written in a familiar vector form: 

d 51 

dt 5uj ' 5uj ds 5f2 

dU 



51 d 51 ^ 51 51 

X uj — —^7^ — Six ^TTT + p X — 

5p 



5n 



s, s')^ (s, s') +^(s,s') 



dU 



[s,s 



ds', 



1.28) 



In order to close the system, one computes the time derivative of ^(s, s') = A ^(s')A(s): 

e(s,s') = -A-i(s')A(s')A-i(s')A(s) + A-i(s')A(s) 

= -cu(s')e(s,s')+e(s,s')u;(s). (1.29) 



Ellis et al. 



Dynamics of charged molecular strands 



18 



This expression is not quite a commutator because different positions s and s' appear 
in u. However, operating with from the left in equation (1.29) gives a proper 
Lie-algebraic expression for the reconstruction of the relative orientation, 

r'iis, s') = uj{s) - Ad5-i(,,,,)cu(s') . (1.30) 

Formulas (1.27) - (1.29) generalize the results in [13] for a fixed filament from 5*0(2) 
to 5*0(3) rotations. 



2 Motion of exact self-interacting geometric rods 
2.1 Problem set-up 

Suppose each rigid conformation of charges is identical and the k-th electrical charge 
is positioned near a given spatial point r through which the curve of base points of 
the RCCs passes. This curve is parametrized by a variable s which need not be the 
arc length. Rather, we take s G [0, L] to be a parameter spanning a fixed interval]^ 
The spatial reference (undisturbed) state for the k-th charge in a given RCC is the 
sum r(s) + r/^(s). That is, r]i^{s) is a vector of constant length that determines the 
position of the k-th electrical charge relative to the point r(s) along the curve in 
its reference configuration. The r]j^{s) specify the shape of the rigid conformation 
of charges. At time t the position of the k-th charge in the rigid conformation 
anchored at spatial position r{s) along the curve parametrized by s may rotate to a 
new position corresponding to the orientation A{s,t) in the expression 

Ck{s) = r{s) + A{s,t)r]f,{s) , where A(s,0) = ld. (2.1) 

This rigid conformational rotation is illustrated in Figure [l} In Mezic's case [13], the 
rotation is in the plane, so that A G 50(2), and there is only one charge, so = 1. 



2.2 Convected representation of nonlocal potential energy 



One part of the potential energy of interaction between rigid conformations of charges 
at spatial coordinates r(s) and r{s') along the filament depends only on the magnitude 
|cm(s') — Cfc(s)| of the vector from charge k at spatial position Cfc(s) to charge m at 
spatial position Cm{s'). This is the Euclidean spatial distance 



(2.2) 



between the k-th and m-th charges in the two conformations whose base points are 
at r{s) and r{s'), respectively. In this notation, the potential energy is given by 



E = E,,Jn,T)- 



dr 




dr 









dsds' 



(2.3) 



^Note: limiting its parametrization to a fixed interval does not mean that the filament is inex- 
tensible. 
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for an appropriate physical choice of the interparticle interaction potential U{d 



k,m I 



and the quantities fi, T (and uj, 7 and p below) are defined in (1.1). The part 
Eioci^, r) represents the purely elastic part of the potential, and is usually taken to 
be a quadratic function of the deformations (f2, F), but more complex expressions are 
possible as well; we shall not restrict the functional form of that dependence. The 
total Lagrangian I is then written as a sum of local lioc and nonlocal Inp'- 

hoc = K{u}, 7) - Eioci^, r, p) and Inp = -E^p , (2.4) 

where K is the kinetic energy that depends only on the local velocities cj, 7. For the 
sake of generality, here and everywhere else below we shall simply consider the total 
Lagrangian to be a sum of the local part lioc{i^,'y,^,T, p), and the nonlocal part 



given by (2.3 ): 



l = lioc{u;,j,n,T,p) + lnp. (2.5) 



The scalar distance dk,m in (2.2) and (2.3) may also be expressed in terms of vectors 



seen from the frame of orientation of the rigid body at a spatial point r(s) along the 
filament, as 

dk,m{s,s') = \Cra{s') - Cfc(s)| 

= \A-\s) (c„(s')-Cfc(s))| 

= \A'\s) {r{s') - vis)) + A-\s)Ais')r,^is') - r,,is)) \ 

=: \K{s,s')+as,s)rt^{s')-r]k{s)\ , (2.6) 

where we have defined the quantities 

k{s,s') := A-\s){r{s') -r{s)) eR"^ and ^{s, s') := A-\s)A{s') e S0{3) . (2.7) 

The first of these quantities is the spatial vector from r(s) to r{s'), as seen from the 
orientation A(s) of the rigid charge conformation located at coordinate label s along 
the filament. The second is the relative orientation of the rigid charge conformations 
located at coordinate labels s and s'. For later use, we record the transposition 
identities, 



which follow from the definition of ^(s, s') in (2.7). 



Remark 2.1 (Left 5*0(3) invariance) Both the body separation vector k(s, s') and 



the relative orientation ^(s,s') defined in (2.7) are invariant under rotations of the 



spatial coordinate system obtained by the left action 

{r{s') -r{s)) ^ 0{r{s') ~r{s)) and A ^ OA , 



by any element O of the rotation group 50(3). 
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Proposition 2.2 (Left SE{3) invariance) The quantities k) G S0{3) x de- 
fined in (2.7) are invariant under all transformations of the special Euclidean group 
SE{3) acting on the left. 

Proof. As a set, the special Euclidean group SE{3) is the Cartesian product SE{3) = 
5*0(3) X whose elements are denoted as (A,r). Its group multiplication is given, 
e.g., in [21] by the semidirect-product action, 

(Al, ri)(A2, rs) = (A1A2, ri + Aira) , (2.9) 

where the action of A G 5*0(3) on r G is denoted as the concatenation Ar and 
the other notation is standard. For the choice 

(Ai,rO = (A,r)-i(s) and (A^, r^) = (A, r)(s') , 

the SE{3) multiplication rule ( 2.9[ ) yields the quantities (^, n) G 5*0(3) x M? as 

(A,r)-i(s)(A,r)(s') = (^(., s'), .c(s, .')) ■ (2.10) 

This expression is invariant under the left action (A, r) (O, i')(A, r) of any element 
(O, v) of the special Euchdean group SE{3). ■ 

Remark 2.3 The 5*ii^(3) setting will be especially important to the development of 
the Lagrange-Poincare formulation of the dynamical filament equations in Section [8j 

Next, let us define the following 5*i^^(3)-invariant quantities, where prime denotes the 
derivative with respect to s and dot is the derivative with respect to t: 



(2.11) 





= A- 


'A' 


G50(3), 


OJ 


= A- 


-lA 


Gso(3), 


r 


= A- 


-v 


G M^ 


7 


= A- 


-If- 


G M^ 


p 


= A 


-\r - ro) 


G 



Hereafter, we shall choose ro = to recover the bundle coordinates (1.1). 

Remark 2.4 Note that here A, r, Q, u, F, 7, p are interpreted as functions of the two 
variables s and t. It will be important to see these variables as time-dependent curves 
with values in functions spaces. For example, we can interpret A[s,t) as a function 
of space and time 

(s, t) G [0, L] X M ^ A(s, t) G S0{3), 
or we can see A as a curve in an infinite dimensional Lie group 

t G M ^ A(-, t) G J^{[0, L], S0{3)), 
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where J-'{[0, L], S0{3)) denotes the group of smooth functions defined on [0,L] with 
values in SO (3). 

This observation is fundamental and leads to two different geometric approaches to 
the same equations: the affine Euler-Poincare and the covariant Lagrange- 
Poincare approaches. 



Remark 2.5 Since A G 50(3), one finds that 



dr 













Wis) 



and the nonlocal part of the potential energy in (2.3) reduces to 

Enp = -Y.\ [[ U{dk,m{s,s'))\r{s)\\r{s')\dsds'. 



(2.12) 



(2.13) 



Remark 2.6 Everywhere in this paper, we shall assume that the nonlocal part of 
the Lagrangian /„p is a function or functional of F, ^ and n. It could, for example, 
be expressed in the integral form 

/„p(e,Ac,r) = jj f/(e(s,s'),K(s,s'),r(s),r(s'))dsds' (2.i4) 



or be a more general functional. In this work, we shall consistently use formula (2.14) 
to make our computations more explicit, although of course the methods would apply 
to more general functionals. Clearly, expression ( |2.3 ) is a reduction of (2.14) obtained 
when the energy of the system of charges is a (half)-sum of interactions between all 
charges. This happens, for example, when investigating electrostatic or screened 
electrostatic charges in a linear media. 

Even though the expression Inp = /„p(^, k, F) is rather general, it is interesting to 
note that physical systems exist whose nonlocal interactions do not satisfy that law. 
For example, the electrostatic potential around a DNA molecule immersed in a fluid 
satisfies the nonlinear Poisson-Boltzmann equation and finding the potential in that 
case is a well-known problem for supercomputers [22] . If we could somehow explicitly 
solve this equation - which is impossible - we would be able to write a more general 
Lagrangian. In general, to apply our theories to this problem we would have to couple 
our methods to a numerical solution of the Poisson-Boltzmann equation at each time 
step. We shall also note that for the case of linearized Poisson-Boltzmann equation, 
we can solve the equation exactly in terms of the screened electrostatic potential 
U{r) = e~^'^/r and the expression (2.14) holds. 



2.3 Kinematics 

We first define auxiliary kinematic equations that hold without any reference to dy- 
namics. We call these advection relations, in order to distinguish them from the 
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dynamical equations (derived later) that balance the forces determined from the 
physics of the problem. In contrast, the advection relations hold for all strands, 
irrespective of their dynamic properties. 

In order to derive the first set of advection relations, we compute the time and space 
derivatives of p{s,t) = A~V(s,t). First, the s-derivative along the filament is given 
by: 



and hence equations (2.11) imply 



p' = -Qp + r = -i7xp + r. 

Next, the time derivative is written as, 

o= -A-1AA-V + A-V, 



and equations (2.11) yield the formula, 

p = -up + 7 



-LJ X p + 7 . 



(2.15) 



(2.16) 



(2.17) 



The next set of advection relations is derived by the equality of cross-derivatives with 
respect to t and s for any sufficiently smooth quantity. First, we use the fact that 
dgdtr = dtdgV. Equality of these cross-derivatives implies the relations. 



7 



-f2 X 7 + A^^f' 



and 



r = -u; X T + A-^r' . 
The difference of the last two equations yields the following relation 

r + a;xr = 7' + f2x7. 



(2.18) 



As we shall see later, the latter is a type of zero- curvature relation. Similarly, 
equality of cross-derivatives dgdtA = dtdgK yields the other advection relation. 



S7 — Lj' = LJ X f2 . 



(2.19) 



2.4 Remark on the n-dimensional generalization and the use 
of other groups 

The previous setting may be generalized to n dimensions and to arbitrary Lie groups. 
This is not only useful for the generalization of charged strands to membranes and, 
more generally, to deformable media; it also gives a more transparent vision of the 
underlying geometric structure underlying the phenomena. 
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Consider the semidirect product 0(^E of a Lie group O with a left representation 
space E. The variables r and A defined above are now functions defined on a space- 
time D X M, where D is a n-dimensional manifold: 

A : G P X M ^ A(s,t) G O, and r : (s, t) G I? x M ^ r(s, t) G 

We will avoid using boldface notation as the functions we consider may be more 
general geometric quantities, not only vectors. As before, "dot" ( ' ) over a quantity 
denotes its time derivative. The derivative with respect to a variable in T) is denoted 
by d; for V = [0,L] this was previously denoted by "prime" ('). The definitions 



(2.11) become 

n = A-^dA : TD ^ , 
u = A-^A -.V 0, 

T = A~^dr -.TV ^ E, (2.20) 



A-V -.V ^ E 



7 

p = A"V -.V ^ E. 

Thus, if we interpret (A,r) as a curve in the group J-'(T>, 0(§)E), the previous defi- 
nition can be rewritten as 

K7) = (A,r)-i(A,r), 
{n,T,p) = c{{A,r)-'), 

where c is defined by 

c(A, r) = ((A, r) d(A, r)-\ -r) . (2.21) 

Remarkably, c is a group one-cocycle. Thus, the previous definition simply says 
that {Q,r,p) are affine advected quantities with zero initial values. This observation 
strongly suggests a relation with the affine Euler-Poincare theory developed in the 
context of complex fiuids in [17j . 

On the other hand, if we interpret (A, r) as a section of the trivial principal bundle 

{VxR)x 0(§)E ^VxR 

over spacetime, definition ( |2.20[ ) simply says that the variables (fi, to, p) are obtained 
by reduction by the subgroup O of the first jet extension of (A,r). This, in turn, 
leads to a relation with the covariant Lagrange-Poincare reduction for field theories 
developed in ^j. Note that by choosing the one-dimensional interval V = [0,L], 
the Lie group O = SO (3) and left representation space E = M.^, one recovers the 
advection of charged strands discussed earlier. 



Remark 2.7 Generalizing to higher dimensions reveals certain distinct aspects of 
the underlying geometry of the problem that are not distinguished in considering the 



Ellis et al. 



Dynamics of charged molecular strands 



24 



particular case of the charged strands. For example, in the case of charged sheets or 
charged elastic deformed media, X> is a domain in M."^, with n = 2 or 3, respectively, 
so the coordinate s has several dimensions. Then, T should be considered as a set 
of vectors Fi, . . . ,r„. Likewise, for the problem of flexible strands of rigid charge 
conformations the distinct objects E and o both coincide with M^. This coincidence is 
removed in higher dimensions and thereby clarifies the underlying geometric structure 
of the theory. 



3 Derivation of the equations of motion 

In this section we shall derive the convective equations of motion for a charged strand 
from two different, but equivalent, viewpoints. The first derivation is based on the 
classical Hamilton-Pontryagin (HP) approach in control theory (see, for example, 
|24j). The second derivation is based on the Euler-Poincare (EP) approach, modified 
to include additional terms describing nonlocal contributions. We shall present both 
methods in this section. 



The Hamilton-Pontryagin Theorem |3.1| elegantly delivers the key formulas for the 
Euler-Poincare equations and leads efficiently to its Lie-Poisson Hamiltonian formu- 
lation. Perhaps surprisingly, the HP theorem produces these results without invoking 
any properties of how the invariance group of the Lagrangian acts on the configura- 



tion space (a manifold) and leads directly to the equations of motion (3.3) and (3.4). 
The equivalent alternative EP derivation of these formulas does explicitly involve the 
action of the Lie group on the configuration space and is, therefore, slightly more 
elaborate than the HP theorem. This elaboration invokes the Lie group action on 
the configuration space and thereby provides additional information. In particular, 
the EP approach reveals how the Lie group action on the configuration space induces 



the affine structure of the EP equations (3.29) and (3.32). The alternative EP ap- 



proach also yields information that explains precisely how the canonical phase space 
(the cotangent bundle of the configuration manifold) maps to the Lie-Poisson space 
associated to the action, which is the dual of the Lie algebra of symmetries via the 
momentum map defined by the infinitesimal affine Lie algebra action. We explore in 
detail the EP route in this paper because it explicitly reveals the role of the Lie group 
action in symmetry reduction. In Section [6] it will be shown that the derivation of 
the EP equations and of the associated variational principle are corollaries of general 
theorems for systems whose configuration space is a Lie group. The complementary, 
but less transparent, HP route reveals other perspectives and results whose abstract 
general formulation will be explored in future work. 
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3.1 A modified Hamilton-Pontryagin approach 
3.1.1 Filament dynamics 

We begin with tlie Hamilton-Pontryagin approach applied to the case when the La- 
grangian includes only the local part, so I = /zoc(<^, 7, f^, T, p). In order to simplify 
the formulas and avoid extra factors in the integrals, we shall implicitly incorporate 
the dependence of the nonlocal potential on F 



A-^r'. See ([3^ and (|332| be- 



low for the explicit formulas for the pairwise potential, that gives multiplication by 
|F(s)||F(s')| in the integrals. 

Inspired by the classical Hamilton-Pontryagin approach, we introduce Lagrange mul- 



tipliers for the holonomic constraints that impose the defining relations (2.11) for the 
five quantities (cj, 7, fi, F, p). 

Theorem 3.1 (Hamilton-Pontryagin theorem for filament dynamics) 

The equations for filament dynamics arise from the variational principle 6S = with 
action S given by 

s = j i{u},j,n,T,p)dt + jj (^TT ■ (^A-^A-u^ +n- {A-^A' -n) 

+ R ■ (A- V - p) + p ■ {A-^r - 7) + M • (A- V - F) ^ & dt. 



These equations are 



51 



R = 0, 



51 



77 = 0, 



51 



51 



51 

5p 5uj 5il 5j 5T 

7r + ujX7r + II' + flxII + jxp + TxM + pxR = 0^ 



M = 0, 



and 



/:i + Ljxp + M'-hrixM-R = 0, 

together with the constraints, 

A-^A = uj, A-^A' = Q, A^r = p, A'^r = j , A" V = F. 
We begin by computing the variations of the quantities appearing in the action S. 



Lemma 3.2 The variations of the quantities in A and r of the formulas in (2.11) 
are 

9S 



5 i^A-^A 
5 {A-'A') 



dt 
S 



+ 



A-^A,S , 

a-^a',e1 , 

5 (A-^r) = * - S (A-V) , 
5 {A-H) = * - S {A-^r) + (a^^a) * , 
5 (A-^r') = - S (A- V) + (A-^A') * . 
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Proof. We calculate the variations directly, one by one. First we have, 
6 (a^^a) = -A-^6A (a^^a) + A-^6A 

= -A-^6A (a^^a) + (A-^5A)- + (a^^a) (A'^M) 



~dt 



+ 



A-iA,S 



Similarly, for the variation of A A' we have 



A-^A'^S 



6 {A-^A') = S 

Now we consider the variation of A^^r, which is given by 
S{A-^r) = -A-^6AA-^r + A-^6r 

= - {A-^6A) (A-V) + {A-^r)- + (a^^a) (A-V) 

= * - S {A-^r) + (a^^a) * . 

A similar argument yields the variation of A^^r', 

6 (A-V) = - S (A-V) + (A-^A') * . 
Finally, the variation of A"^r is given by, 

S (A-V) = -A-^5AA-V + A-^Sr = * - S (A'V) 
and all the formulas in the statement are proved. ■ 



We may now use these identities to prove the Hamilton-Pontryagin Theorem 3.1 for 
the equations of filament dynamics. 

Proof. The main results from this Hamilton's principle arise from the following iden- 
tities, written in terms of the skew-symmetric 3x3 matrix S = A^^6A G 5o(3) ~ 
and the vector ^ = A~^6r G M^. Variations with respect to the Lagrange multipliers 
impose the expected defining relations for the five quantities {p, fi, 7, F}. The 
conjugate variations give 



61 ^ SI SI ^ 61 

R = 0, —-77 = 0, — -n = 0, — -^ = 0^ 
oa; oil 07 



61 



6p ' 6u; ' 6n ' 6j ' 6T 

Finally, the variations proportional to S and ^ yield the filament equations 
7r + LJXn + Il' + nxIl + jXfi + TxM + pxR = 0, 

and 

fi + u;xp + M' + nxM-R = 0, 

respectively. ■ 



M = 0. 



Ellis et al. 



Dynamics of charged molecular strands 



27 



Remark 3.1 The Hamilton-Pontryagin approach used here also allows nonholo- 
nomic constraints to be imposed on the motion of the strands, if one desires. See 
|21j for a discussion of nonholonomic constraints using the Hamilton-Pontryagin ap- 
proach. 



3.1.2 Nonlocal potential 



For the nonlocal potential (2.13) we may form a Hamilton-Pontryagin variational 
principal in a similar fashion. In this case, the action Snp is given by 

Snp = J Inpi^, K, T)dt + jj'm- {A-\sy{s) - T) dsdt 

+ III (X- (A-^(s)A(s') -()+K- {A-\s) (r(s') - r(s)) - h) jdsds'dt 



Lemma 3.3 The additional variational formulas needed for calculating the equations 
of motion are given by 

A-\s')Ais) {5 {A-\s)Ais'))) = -AdA-i(,,)A(.)S(s) + S(s') , 
6 {A"\s) {r{s') - r{s))) = -^{s)A~\s) {r{s') - r{s)) 

+ A-i(s)A(s')*(s') + *(s), 

where the independent variations are defined by 

^{s) = A-\s)6r{s) and ^{s) = A-\s)6A{s) . (3.1) 

Proof. The first variational formula is calculated directly, as 

A-\s')A{s) {6 {A-\s)A{s'))) = A-\s')A{s) {A-\s)6A{s')) 

-A-\s')A{s) {A-\s)6A{s)A-\s)A{s')) 

= -AdA-i(y)A(.)S(s) + S(s'). 

The second variational formula follows similarly from a direct calculation, 

6 {A-\s) {r{s') - r{s))) = -A-\s)6A{s)A-\s) {r{s') - r{s)) 

+A-\s) {6r{s') - 6r{s)) 
= -^{s)A-\s) (r(s') - r{s)) + A~\s)A{s')^{s') + 

which proves the lemma. ■ 

Theorem 3.4 The equations that arise from the variational principle with the non- 
local action 




Snp = III U{^, K, T)dsds'dt + JJ ^ - (A"^(s)r'(s) - r) dsdt 

'x ■ {A-\s)A{s') ~^)+K- {A-\s) (r(s') - r(s)) - k) ) dsds'dt 
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are given by: 



dU dU dU 

^=ae' ^=d^' "^^df 



and 



Txm = J {^{s,s')X{s\s)-X{s,s')C\s,s') + K{s,s')xH,{s,s'))ds', 

m' + nxm = J (^(s, s')K{s', s) - K{s, s')) ds' , 
together with the constraints, 

e = A-\s)A{s') , K = A-\s) {r{s') - r{s)) , T = A-\sy{s) . 



Proof. The proof is obtained by substituting the variations given in Lemma [3.3| into 
the Hamilton's principle for the action in the statement of the theorem. Variations 
in X, 1^ and m yield the constraints, 

e = A-\s)A{s') , ^ = A-\s) {r{s') - r{s)) , T = A-\sy{s) . 

Variations in ^, k, and F yield the relationships 

dU dU dU 

X = —— , K = —— , m = -— . 
dC dn' dT 

Finally, the variations proportional to S(s) and ^{s) yield 

Txm = J {^{s,s')X{s',s) - X{s,s')C\s,s') + K{s,s') X K,{s,s')) ds' 

and 

m' + nxm = j (^(s, s')K{s', s) - K{s, s')) ds' , 

respectively. ■ 

We may combine these nonlocal terms and the local part of the equations to produce 
the full set of equations. These are given by 

Tr + Li;x7r + n' + r2xn + 7X X (M + m) + p X R 

= j {K{s,s') X k{s,s') + Z{s,s'))ds' 

and 

fi + u;xn + {M + my + nx{M + m)-R = J (^(s, s')K{s', s) - K{s, s')) ds' , 
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where one defines 

Z{s, s') := ^s, s')X{s', s) - X{s, s')r\s, s') , (3.2) 

denoted as Z since the right hand side of this equation is in so (3). 
We may now use these functional-derivative relations to express the equations of 
motion in terms of the reduced Lagrangian, / = hoc + Inp- The functional-derivative 
relations obtained in the Hamilton-Pontryagin approach are 



R = 


Shoe 
6p ' 


77 = 


Shoe 

du) 


n = 


Shoe 

5n ' 


= 


Shoe 
6 J 


M = 


Shoe 

6T ' 


X = 


Slnp 


K = 


6lnp 

6k, ' 


m = 


5T 



Substituting these relations into the equations of motion above gives the following 
equations of motion for the charged strand. 



6uj 



6n 



Shoe , S {hoe + ^np) 

— X 7+ ^ 

07 



6r 



, ^hoe 

xT + —— X p 
dp 



/ o , N ^hoe , / o , \ ^ {hoc ~(" ^np) 

{dt + ujx) —— + {ds + fix) — — ^ 

07 oi 



Shoe^ [ ,^dU dU 



[s,s')^ds'. (3.4) 



The term Z{s, s') is the contribution from the nonlocal part of the Lagrangian that 
we have sought. 

Remark 3.2 



The dynamical equations (3.3) and (3.4) must be augmented by the advection 



conditions (2.16), (2.18) and (2.19) in order to close the system 



The resulting system of equations describes an elastic filament with two non- 
local additional components, or degrees of freedom, compared to the ordinary 
Kirchhoff filament (to which the system reduces, when p and ^ are absent). 
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• The two additional (nonlocal) degrees of freedom in p and H (with H = 
^~^d^/ds) will produce an important effect that will distinguish the behavior of 
this system from that of the ordinary Kirchhoff filament. Namely, the presence 
of the two additional equations for p and H raises the order of the equation set. 
In turn, the increase in differential order of the system will produce additional 
modes of excitation for the waves that will propagate along the filament when 
the system is linearized around the static solutions. 



Summary. Equations (3.3), (3.4), (2.16), (2.18) and (2.19) represent the generaliza- 
tion of the Kirchhoff model that we have sought. As we shall see in the next section, 
under a certain transformation of variables this model reduces to a conservation law 
formulated in terms of coadjoint action on se(3) Lie algebras. 



3.2 A modified Euler-Poincare approach 



The Euler-Poincare approach is based on applying Hamilton's variational principle to 
the symmetry-reduced Lagrangian and constraining the variations properly. While 
this will be yet another way of deriving equations (3.3), (3.4), we believe that such 
a "bare hands" derivation will benefit understanding, as it represents a direct and 
explicit derivation of those equations of motion. See [IB] and [25 for an introduction 
to the classical Euler-Poincare approach. Some calculations in this section overlap 
with those in Section 3.1, Nonetheless, we have chosen to present them here for 



completeness of exposition. 



3.2.1 Variations: Definitions 

Let us compute variations of p, u:, 7, f2 and F. We proceed by first computing, 

6p = -A"^5AA~ V + A-^6r = -Sp + * = -Sxp + * = pxS + *, (3.5) 

where we have defined the variational quantities 

S = A-i5A, (3.6) 
^ = A-^6r. (3.7) 

Next, we compute the space and time derivatives of S and ^ along the curve. We 
have the space derivative, 

— = -A-^A'A-^5r + A-^5r' = + A-^5r' = -17 x * + A-^5r' , (3.8) 

OS 

and the time derivative, 

— = -A-^AA-^6r + A-^6r = -a;* + A-^6r = -w x * + A'^Sr . (3.9) 
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A-'A'A-'6A + A-'6A' = + A'^dA' 



Analogously, for the space derivative of S, 

ds 

while the time derivative of S is computed as follows: 

cm 

~dt 

Now we are ready to compute the variations ^7, 6r, 600 and 6Q. The first of these is 

9* 



-A-^AA-^6A + A-^6A' = -uT. + A'^dA . 



(3.10) 



(3.11) 



5-f = -A-'5AA-''p+ A-'5p = -S7 + cu* 



dt 



use (3.9) 



so in vector form, 



9* 



(3.12) 



Likewise, 



5T = -A'^6AA-^p' + A-'6p' = -SF + fi* + 



ds 



use (3.8) 



which has the vector form, 



5r = rxs + rix* + 



9* 

ds 



(3.13) 



Next, 



5uj = -A-^6AA-^A+ A-'6A = -Eu + ojE + ^ = [uj,E] + ^ 

dt dt 



use (3.11 ) 



so expressing these formulas in terms of vectors yields 

9S 



(5ci; = a; X S + 



dt 



(3.14) 



Finally, 



dT. dT, 
Sn = -A-^6AA-^A' + A-^5A' = -Sfi + + — = [r]. Si + — , 

ds ds 



use (3.10) 



so, again, expressing in terms of vectors leads to 



6fl = flxj: 



dE 

ds 



Finally, the variation of ^(s, s') is given by 

{r'S^s, s')) = -Ad5-i(,,,,)S(s) + S(s') . 



(3.15) 



(3.16) 
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3.2.2 Derivation of the equations of motion 

Suppose now we want to compute variations of the reduced energy Lagrangian I which 



is a functional of (p, 7, F, cu, From (2.6) we see that 



=\A^\s)r{s,t) - A-\s,t)r{s' ,t) + 'n,,{s) - ^{s, s')'n^{s')\ 

= \k{s, s') + r],{s) - s')r]^{s')\ , (3.17) 



where we have defined 

k{s, s') = A-\s) {r{s') - r{s)) = ^s, s')p{s') - p{s) 

The variation of k, is then given by 

6k{s, s') =S(s)k(s, s') - *(s) + ^(s, s')^{s') 

= S(s) X k{s, s) - *(s) + e(s, s')*(s') . 



(3.18) 



(3.19) 



Let us first define the Lagrangian / as the sum of a 'local' part lioc and a nonlocal 
part Inp, according to 

/(a;,7, n, T, p, ^, k) := + Lp 

= kUu;, j,n,r,p) + jju {k{s, s'), ^s, s'), r{s), r{s')) dsds' . (3.20) 

Note. From now on, we assume that the nonlocal part of the potential energy f/ is a 
function of the two variables k,{s, s') and ^(s, s'), as well as F, since s is not necessarily 
the arc length. In particular, for a potential energy depending on the distance dk^m, 



the variables k and ^ enter in the linear combination defined by (3.17). In principle. 



the potential energy could have chosen to be an arbitrary functional of A~^(s)r(s), 
A~^(s)r(s') and ^(s, s'). Euler-Poincare methods would be directly applicable to these 
functionals as well. 

The equations of motion are computed from the stationary action principle SS = 0, 



with S = f Idt and / = //oc + ^np in equation (3.20). We have 



5S 



61 



loc 



+ 



6p 

Shoe 

In 



6p 

6Q 



61 



loc 



5j 

6lnp 

6k, ' 



, ^7 



6{lloc + I 



np J 



ST 



6r 



Shoe 
5uj 



5uj 



(3.21) 



5k 



where (■,■) = /( ■ , ■)Rids represents pairing in the filament variable s. We may 
now substitute 5p from (3.5), ^7 from (3.12) and 5Q. from (3.15). We have 



'51 



loc 



6p 



6p 



51 



loc 



5p 



p X S + * 



51 



loc 



5p 



X p, S 



51 



loc 



5p 



(3.22) 
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For 5k, we obtain 



+ 

Next, 



's{ h + Q \^/ s{k + Q r X s + » X ^ + ^ 

oT / \ oT OS 

S{lloc + lnp) v\ ^ / K^loc + ^np) „ S{lioc + lnp) . , . 

and 

^7 ) = ( ^ — ,7xS + uJX* + 



^7 ' / \ ^7 ' dt 



X 7 , S ) + ( X a; - ' * ) • (3-25) 



^7 ' / \ ^7 S-y 
Variations in uj and Q give, respectively, after integrating by parts. 



and 



Shoe /Shoe ^ ^ (9S \ / ^ <9 Slloc ^\ 

Finally, one computes the variations in ^ as follows: 

= / (r'(^,5')^(5,50,-Ad^-i(.,.oS(^) + S(s')\ ds', (3.28) 

\ / so(3) 

where (-, ■)3o(3) : so(3)* xso(3) — > M is the real-valued pairing between the Lie algebra 
50(3) and its dual 50 (3)*. 
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Substitution of (3.22), (3. 25), and (3.26) gives an expression for 6S that is linear in 



I] and ^. Collecting those terms when imposing 5S = implies from the term 
proportional to S that: 



d 6lioc , Sli 

oc 

\- LJ X 

00l> 



dt Suj 



d SI 



loc 



ds 5n 



n X 



Shoc\ Slioc _ . S{lioc + lnp) 



6n 



^7 



X 7 + 



ST 



, Slloc 
+ X P 



dU 



Sp J \di^ 

where the term Z{s, s') is the vector given by 



s,s') X k{s,s') + Z{s,s') ds' , 



X r 



(3.29) 



T 



dU 



s,s')e{s,s') 



(3.30) 



which is the same quantity that we found using the Hamilton-Pontryagin approach. 



Formula (3.30) is computed from the variation in (3.28) as follows 



dU 

r\s, s') — {s, s') , -Ad5-i(,,,,)S(s) + S(s') ) dsds' 



so(3) 



-Ad^-i(,^,/)^^(s, s')^{s, s') + ^^(s', s)^{s', s) , S(s) ) dsds' 



so(3) 
T 



II /-as,s')eis,s')^is,s')e{s,s')+as,s') , S(s)\ dsds' 

\ / so(3) 

^^(s,s')eis,s')+as,s')(^is,s')\ ,S(s)\ dsds'. (3.31) 

Here, we have used the fact that ^"^(s, s') = C,~^{s, s'), and ^(s', s) = ^~^{s', s). 
Next, we collect the terms proportional to ^ in order to close the system. We find 

d Slloc , Slloc\ , / d S{lloc + lnp) , ^ S{lloc + lnp) 
(Jj X — 1 + ^ ; — + f2 X ^ — 



dt (57 57 / V ds ST 

Shoe 



ST 



Sp 



Remark 3.3 Equations (3.29) and (3.32) obtained by the Euler-Poincare approach 
recover equations (3.3) and (3.4), respectively, from the Hamilton-Pontryagin ap- 
proach. 
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4 Conservation laws 

In order to elucidate the physical meaning of the somewhat complex-looking equations 
(3.3) and (3.4), we shall write them explicitly as conservation laws. For this purpose, 
we invoke the following identities valid for any Lie group G. Given a smooth curve 
g{t) E G, rj E g, and /i G g*, we have 



Ad3-i(t) ■^Adg(^t)V = ad^(t)^ , (4-1) 
d 

Ad;w^Ad;-i(,)/i = -ad:(,)/i , (4.2) 

where a(t) = g~^g(t) G and Ad* denotes the coadjoint action of G on q* defined by 
(Ad*/i, T]) := (/i , AdgT]). Formula (4.2) generalizes to a curve fi(t) as 



d 



(4.3) 



To derive the conservation form of equations (3.3) and (3.4) we need to consider the 
group G = SE{3) whose elements are denoted hj g = (A,r). Consider the function 
{A{s,t),r{s,t)) defined on spacetime. Then we have 



a 



(A,r)-i(A,r) = (A-iA,A-ir) = (a;,7). 
Recall that the infinitesimal coadjoint action on se(3)* is 



ad 



(u; X + 7 X /3, a; X /3) . 



(4.4) 



(4.5) 



Then, using equations (4.3) and (4.5) for the temporal dual Lie algebra elements 
(yti, (3) = {61/6LJ , 51 /5j) yields 



Ad 



d_ 

(A,r)^ 

d_ 
Ft 



Ad 



(A,r)-i 



loc 



Shoe Shoe 
5uj ' ^7 



Hoc 

5uj ' ^7 

, , 5lloc , 5lloc Slloc 

+ \ (jj X — h 7 X — — , u; X 



(4.6) 



For the derivative with respect to curve parametrization s, we need to remember that 
the nonlocal part of the potential depends on T as well. Thus, we have 



Ad 



d_ 

d_ 
ds 



Ad(A,r)- 

5lioc 51 



5lloc 5{lloc + lnp) 



6n 



loc 



6n ' 6r 



6r 



6n 



6r 



fl X 



5ho 



6r 



(4.7) 



Some additional identities derived below will be needed in treating the nonlocal part 
of the potential. 
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First we deal with the nonlocal term by referring to equation (3.2). This can be 



expressed as a formal derivative of the nonlocal part of the potential with respect to 
Lie algebra elements f2 and T as follows. Note that there are only two free variations 
I] = A~^6A and ^ = A~^6r. On the other hand, the nonlocal part of the Lagrangian 
depends on three variables p,^, and T. Thus, there must be a relation between the 
partial derivatives of the nonlocal part of the Lagrangian and the total derivatives 
with respect to F and fl. This relation is computed as follows. 

Upon identifying coefficients of the free variations Sx = A~^6A and ^ = A^^Sr, the 
following identity relates different variational derivatives of the nonlocal potential Inp' 



— 1 ^^np 



61 



np 



6k, 



6k) + 



61 



np 



6r 



61 



np 



6T 



6T 



Tot 



6T 



61 



np 



6^ 



6Q, 



(4.8) 



We will discuss this point in detail in ^6.3.2, Here, the subscript on (OIto* denotes 



the total derivative with respect to F. Using expressions (3.16) for ^ ^5^, (3.19) for 



6k, (3.15) for 6n and (3.12) for 5F, then collecting terms proportional to the free 



variation S yields the following identity, which implicitly defines 6lnp/ 6Q, in terms of 
known quantities. 



d 61 



np 



ds 6n 



6lnp 



dU . 
Ok 



6n 

[S, s') X K 



(4.9) 



s,s') ds' + j Z(s,s') ds'. 



where we have defined Z(s,s') according to (3.2). Likewise, identifying terms multi- 
plying ^ gives 



d 61 



np 



ds 6T 



n X '^^"^ 



Tot 



6T 



d 61 



np 



Tot 



ds 6T 



+ 



dU , 



6ljip 
ds'. 



(4.10) 



Therefore, we conclude that equations (3.3), (3.4) are equivalent to the following 



equations expressed on 5e*(3) in conservative form using variations of the total La- 
grangian, / := lioc + Inp- 



d_ 
di 



[6^' 6^) 



d_ 
ds 



61 61 



\6n' 6T 



Tot / J 



61 



61 



Here, the components of 



AW \jp''P^6p 



61 61 



(4.11) 
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represent, respectively, the spatial angular momentum density and the spatial linear 
momentum density of the strand, whose center of mass lies along its centerline. The 
components of 



are the external torques and forces. (See (1.23) for the last simplification.) As men- 



tioned above, only external forces arising from potentials are considered in this paper. 
In principle, more general non-conservative forces and torques can be considered as 
well, but we shall leave this question for further studies. 

Remark 4.1 For future reference, it is advantageous to write out the conservation 



law (4.11) in convective form as 



, 51 ^ , 51 51 ^ 51 51 

9t + LJx — + 9, + f2x — + px — + rx — +7X — = 0, 

5. 5n 5p 5T 5j ^,.,2) 

idt + u;x) — + {d, + nx) — - — =0. 
07 oi op 

Here we have defined the total Lagrangian / := lioc + Inp, and all the variational 
derivatives are assumed to be the total derivatives. Note that in these equations 
coincide precisely with the equations for the purely elastic filaments derived in [T2] . 



We note that the variations with respect to fl and F are computed implicitly in (4.9 



4.10). To actually use these equations to explicitly describe nonlocal interactions. 



we must expand the derivatives with respect to ^ and k in (4.12). However, we 
emphasize again that it is interesting that nonlocal interactions can be expressed so 
as to formally coincide with the equations for the purely elastic motion. See §6.3.2 
for a detailed discussion of this point. 



5 Hamiltonian structure of the strand equations 

It is useful to transform the Lagrangian dynamical equations into the Hamiltonian 
description, both to relate these equations to previous work on elastic rods and to 
elucidate further their mathematical structure. We start by Legendre transforming 
the total Lagrangian / to the Hamiltonian, 

h{p,f3,n,r,p)= [{p-u; + f3--f)ds-liu;n,^,r,p), (5.1) 



where u: , 7 are determined from the relations p = 51/ 5uj and (3 = 51/5^ upon assum- 
ing that / is hyperregular. Then, equations ( 2.17[ ), (2.18), (2.19), and ( |4.12 ) may be 
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expressed in Lie-Poisson form with three cocycles as 







fix 


f3x 


{ds + nx) 


rx 


px 






I3x 








{ds + nx) 


-Id 


n 


















r 




Tx 


(ds + nx) 











p 




px 


Id 














5h/5n 




5h/5(3 




5h/5n 




6h/6T 




5h/5p 



d_ 
dt 



(5.2) 

Note that a; = 6h/6p and 7 = 6h/6f3. The affine terms dg and Id arise from the 
cocycle appearing in the definition of the variables fl,r,p in (1.1); see also (2.20). 
These equations produce the affine terms located in the matrix elements {p, f2}, 
{/3,r}, and {/3, p}. 

This Hamiltonian matrix defines an affine Lie-Poisson bracket on the dual of the 
semidirect product Lie algebra 

J^(/,se(3))(s)J^(/,se(3) x R^), 

where se(3) =so{3)(i)R\ I = [0,L], and 

{p, (3) e se(3))* and (12, T, p) E ^{1, se(3) x 

The associated affine Lie-Poisson bracket reads 



{f,g}Xp,f3,ft,r,p) 



- / n 



- / n 



^x^ 
6p 6p 

^x^ 
6fl 6p 

5T 5(3 

^x^ 
5T 6p 

6p 6p 



/3- 



5/ Sg_ 
5(3 ^ 5p 



Sg_ SI 
6(3 ^ 6p 



60, Sp 

5T 6(3 

^x^ 
6T 6p 

6p 6p 



(5.3) 



+ 



5/ 



"5p 5T '5(3 6p 6(3 



59_ 
6(3 
Sf 



a , ^9 6f 6g 

— ■ 8^ ■ 8^ ■ — . 

6n 6p 6T "6(3 6p 6(3 

The first line represents the Lie-Poisson bracket on the Lie algebra JF(J, se(3)). The 
first five lines represent the Lie-Poisson bracket on the semidirect product Lie algebra 

J^(/,se(3))(s)J^(/,se(3) x R?). 

The last two lines represent the affine terms due to the presence of a cocycle, as 
well as the canonical Poisson bracket in (p, (3). The Poisson bracket (5.3) is an 
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extension to include p of the Poisson bracket for the exact geometric rod theory of 
[12] in the convective representation. Remarkably, from a geometric point of view, 
this Hamiltonian structure is identical to that of complex fluids [T7112S]. The reason 
for this will be explained in detail in Section [6j 



6 The affine Euler-Poincare and Lie-Poisson ap- 
proaches 

This section explains how the equations of the charged strand may be obtained by 
affine Euler-Poincare and affine Lie-Poisson reduction. This proves that the charged 
strand admits the same geometrical description as the complex fluids and spin sys- 
tems. 

We begin by recalling from [T7j the theory of affine Euler-Poincare and Lie- 
Poisson reduction. In contrast to [TT], however, we consider here Lagrangians and 
Hamiltonians that are /e/it-invariant, rather than right-invariant. 



6.1 Notations for semidirect products 

Let be a vector space and assume that the Lie group G acts on the left by linear 
maps (and hence G also acts on the left on the dual space V*). As a set, the semidirect 
product S = G(§)V is the Cartesian product S = G xV whose group multiplication 
is given by 

i9i,vi)i92,V2) = igig2,vi + giV2), 

where the action of g & G on v ^ V is denoted simply as gv. The Lie algebra of S is 
the semidirect product Lie algebra, s = q(§)V, whose bracket has the expression 

ad(5,,„,)(6,f2) = [(6, ^^i), (6,^2)] = (['^i,6],6^2 
where denotes the induced action of q on V, that is. 



d 
dt 



exp{t^)v E V. 



t=o 



From the expression for the Lie bracket, it follows that for (^, v) E s and (/x, a) G s* 
we have 



ad 



(/i, a) = (ad^ fi — V o a, —^a) 



where E V* and v o a E Q* are given by 



d_ 
dt 



exp(t^)a and {v o a,^)g := —{^a,v) 



V, 



t=0 
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and where (•, •)g : x ^ R and (•, •)y :V* xV are the duahty pairings. The 
coadjoint action of S on 5* has the expression 

Ad(3^^)-i in, a) = (Ad*-i fx + voga, go) . (6.1) 

Suppose we are given a left representation of G on the vector space V*. We can form 
an affine left representation 9g{a) := ga + c{g), where c G J-'{G, V*) is a left group 
one-cocycle, that is, it verifies the property 

c{gh) = c{g) + gc{h) , (6.2) 

for all g,h e G. Note that 



d 
dt 



6'exp(to(a) = C« + dc(0 



t=0 



and 

(^a + dc(0,'y)y = (dc^(T;) -i;^a, 
where dc : Q ^ V* is defined by dc(^) := T'ec(^), and dc^ : F ^ 0* is defined by 

6.2 AfRne Lagrangian and Hamiltonian semidirect product 
theory 

Concerning the Lagrangian side, the general setup is the following. 

• Assume that we have a function L : TG x V* ^ M. which is left G-invariant 
under the affine action {vh, a) ^ {gvh, dgia)) — {gvh, ga + c{g)). 

• In particular, if oq G V*, define the Lagrangian L^g : TG — > R by LaQ{vg) :— 
L(vg, Oq). Then is left invariant under the lift to TG of the left action of 

on G, where G"^^ is the isotropy group of Oq with respect to the affine action 9. 

• Define I : q x V* ^ M hy I :— L\gxv- Left G-invariance of L yields 

l{g''^Vg,eg-i{a)) = L{vg,a) 
for all geG,Vge TgG, a e V*. 

• For a curve g{t) e G, let ^{t) :— g{t)~^g{t) and define the curve a{t) as the 
unique solution of the following affine differential equation with time dependent 
coefficients 

d = -^a - dc(0, 

with initial condition a(0) = oq. The solution can be written as a{t) ~ 
Og{t)-^{ao)- 
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Theorem 6.1 In the preceding notation, the following are equivalent: 
i With ao held fixed, Hamilton's variational principle 



6 [\ao{9,9)dt = 0, (6.3) 
J to 



holds, for variations Sg(t) of g{t) vanishing at the endpoints. 
ii g{t) satisfies the Euler- Lagrange equations for La^ on G. 
iii The constrained variational principle 



6 l{^,a)dt = 0, (6.4) 

J to 

holds on Q X V* , upon using variations of the form 

dfi 

5^ = — + [^,7]], 6a = -ria- dc{r]), 

where ri{t) G g vanishes at the endpoints. 
iv The affine Euler-Poincare equations hold on q x V* : 

— — = adf — + — o a - dc . (6.5) 

dt5^ ^5^ 5a \5aJ ^ ^ 

See [T7j for the proof and apphcations to spin systems and complex fluids. Concerning 
the Hamihonian side, the setup is the fohowing. 

• Assume that we have a function H : T*G x V* ^ which is left invariant 
under the affine action {ah, a) ^— {gah,Og{a)). 

• In particular, if oq G V*, define the Hamihonian if^o • T*G — > M by HaQ^Oig) : = 
Hi ag, Oq). Then is left invariant under the lift to T*G of the left action of 
Gl, on G. 

• Define h : g* x V* Rhy h := Hl^^y^v* ■ Left G-invariance of H yields 

h{g-^ag,9g-i{a)) = H{ag,a). 
for aU g e G, ag e T*G, a eV*. 
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Note that the G-action on T*G x V* is induced by the S'-action on T*S given by 

(u, «)) := (g^h, V + gu, ga + c{g)) . (6.6) 

The affine action \l/ appears as a modification of the cotangent hft of left translation 
on S by an affine term. Thus, we can think of the Hamiltonian H : T*G x V* ^ M. 
as being the Poisson reduction of a S*- invariant Hamiltonian H : T*S hj the 

normal subgroup {e} x V since (T*S)/{{e} x V) = T*G x V*. Note also that every 
Hamiltonian H = H{ah, {u,a)), defined on T*S and left invariant under the affine 
action does not depend on the variable u E V. 

Theorem 6.2 Let a{t) G T*^^-^G he a solution of Hamilton's equations associated to 
with initial condition jiQ e T*G. Then {fi(t),a(t)) := {g(t)^^a(t),9g(^t)-i{ao)) G 
Q* X V* is a solution of the affine Lie-Poisson equations on s* : 



d_ 
di 



(fx, a) 



S/j, 



Sh 

5a 



oa + dc 



f5h 



5h 



dc 



5h\ 
5fiJ 



with initial conditions (/i(0),a(0)) = (yUo,cto) 
affine Lie-Poisson bracket on the dual s* 



The associated Poisson bracket is the 



+ ( dc 



SI s_i 
Sfi' Sfi 

s_i 

Sfi 



Sa 



a 

' 5/i Sa 
-/dc 



_sis_i 

5/i Sa 

s_i\ s_i 

SfiJ ' Sa 



(6.7) 



Conversely, given fiQ G T*G, the solution a{t) of the Hamiltonian system associated to 
is reconstructed from the solution {fi(t) , a{t)) of the affine Lie-Poisson equations 
with initial conditions (/i(0),a(0)) = (/io,ao) by setting a(t) = g(t)fi{t), where g(t) is 
the unique solution of the differential equation g{t) = g{t)-^^ with initial condition 
g{Q) = e. 

Proof. See [T7j for the proof and some applications. ■ 



Momentum maps We now comment on the momentum maps at each stage of 
the reduction process. In [T7j it is shown that the momentum map associated to the 



affine action (6.6) is given by 



3:T*S^s\ J(a3, (m,6)) 



{ugg ^ + uob — dc'^{u),b). 



(6.8) 



The proof of this formula uses the general formula for the momentum map on a 
magnetic cotangent bundle with respect to the cotangent-lifted action. In order to 
apply this formula, an adequate fiber translation on T*S used. This fiber translation 



turns out to be equivariant with respect to the affine action and the action (6.6) 
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on T*S as well as symplectic with respect to the canonical symplectic form and a 
magnetic symplectic form on T*S. The above formula for J is then obtained by 
pulling back the magnetic momentum map via the fiber translation. 
One observes that the conservation of J implies that the motion takes place on affine 
coadjoint orbits. 

Note that the Poisson action of G on T*G' x V* does not admit a momentum map 
because the leaves T*G x {b}, b G V*, are not invariant under this action. Given 
ao ^ V*, the momentum map on T*G corresponding to the cotangent lifted left action 
of the isotropy group G"^ is given by 



where 0^ = G g | + dc(^) = 0} is the Lie algebra of G^^. 
6.3 Affine reduction at fixed parameter 

As we will see, the affine reduction theorems recalled above do not apply directly 
to the molecular strand. This is because the Lagrangian of the molecular strand is 
only given for the particular value ao = of the parameter and we do not have a 
concrete expression for La^ when ao 7^ is an arbitrary element of V*. Extending Lq 
by G-invariance only yields a Lagrangian on TG x Oq, where Oq C V* is the orbit of 
the affine G-action on V*. Fortunately, the Lagrangian Lq for the molecular strand 
is invariant under the isotropy group Gq = {g E G \ c{g) = 0} and this turns out to 
be enough for the extension of the affine semidirect product reduction theorem. 

6.3.1 Lagrangian approach 

We consider here the case of a GJ^^ -invariant Lagrangian La^ : TG M for a fixed 
ao G V*, but we do not suppose that this Lagrangian comes from a G-invariant 
function L : TG x l^* — > M. In particular, we do not know the expression of La when 
a 7^ is an arbitrary element of V*. To La^ we associate the reduced Lagrangian / 
defined on the submanifold 



given by /(^, 6'g(ao)) = La^^g ^^). The tangent space at a to O^^ is given by 



Jao : T*G ^ (g^J*, Jao(ag) 




gxOl^Cgx V*, := {^,(ao) | g G G} 



T O' 



{dc{r]) + T]a \ T] e g}. 



(6.9) 



The analogue of Theorem 6.1 in this case is given below. 



Theorem 6.3 Let ao be a fixed element in V* and g{t) be a curve in G with g{0) = e. 
Define the curves ^(t) = g(t)~^g(t) G Q and a(t) := Og^t^iao G V* . Then the following 
are equivalent. 
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i With oo held fixed, Hamilton's variational principle 



S L,,{g,g)dt^O, (6.10) 

J to 



holds, for variations 5g{t) of g{t) vanishing at the endpoints. 
ii g{t) satisfies the Euler- Lagrange equations for La^ on G. 
iii The constrained variational principle 



S [ l{^,a)dt^O, (6.11) 



'to 

holds on Q X O^^ G Q x V* , upon using variations of the form 

df] 

5^^ — + [i,r]], 6a = -r]a - dc{r]) , 

where r]{t) eg vanishes at the endpoints. 

iv Extending I arbitrarily to Q xV* , the affine Euler-Poincare equations hold on 
the submanifold g x O^^ (Z QxV*: 

d 61 ^,61 61 ^Tf^^\ /..ON 

mrr'^^M + 6a''-''' [ra)- ^'-''^ 

Proof. The equivalence of i and ii is true in general. The equivalence of i and iii 
and the equivalence of iii and iv can be shown exactly as in the standard case, that 
is, the case when / is defined on the whole space q x V*. The only minor difference 
occurs when I is differentiated with respect to the second variable. In this case the 
functional derivative Sl/Sa e F is replaced by the tangent map d2l{^, a) e T*Ol^ and 
one observes that 




d2/(^, a)-6a = ( — , 6a ) , for all 6a e TaO[ 



ao 



for any extension / of / to g x \/*. Note that 6a = —rja — dc{ri) G TaO^^ foi rj E rj and 
that any vector in T^C^^ is of this form. From now on we denote also by /, instead 
of Z, an arbitrary extension of Z. ■ 

Remark 6.4 (The case ao = and the charged strand) For the charged molec- 
ular strand we will need to choose Oq = 0. In this case the isotropy group is 
Gq = {g e G \ c{g) = 0}. Given a GQ-invariant Lagrangian Lq : TG — > R, the 
reduced Lagrangian I is defined on qx Oq by 

l{Lc{g-')) = L,{gi). 
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It will be sufficient to restrict to Lagrangians for simple mechanical systems with 
symmetry, that is, of the form -Lo(fg) = K{vg) — P{g), where K is the kinetic energy 
associated to a Gg-invariant Riemannian metric on G and the potential P is Gq- 
invariant. In this case, the reduced Lagrangian is 

l{^,c{g~')) = K{gO-P{9). 

Note that the right hand side of this expression is well defined on g x Oq, that is, it 
depends on g only through c{g~^). Indeed, c{g"^) = c{h~^) if and only if 6g-i{0) = 
^h-i(O), which means that hg^^ G Gq. Therefore, P{h) = P{{hg~^)g) = P{g) by left 
GQ-invariance of P. For the kinetic energy the same argument works since the metric 
is Go-invariant. 

Thus we can write Lo^Vg) = K{vg) — E{c{g^^)) for the function E : V* ^ M. uniquely 
determined by the relation P{g) = E{c{g^^)). In this case, we have 

li^,c{g-')) = KigO-Eicig-')). 

For the Lagrangian of the charged molecular strand the potential energy is the sum 
of two terms, one of which, denoted by Eioc, explicitly depends only on c{g~^) and the 
other, denoted by Enp, does not have a concrete expression only in terms of c{g~^) but 
it is Gg-invariant. In addition, for the charged molecular strand the kinetic energy 
metric is not just Gg-invariant but G-invariant which then implies that it is only a 
function of G 0. For the molecular strand the Lagrangian is of the form 

Loivg) = K{vg) - EloMg-')) - Er^d), cig~')), 

where C is a Gg-invariant function defined on G and the reduced Lagrangian is 

K^, cig-')) = KiO - EioMg-')) -Er^Ma), c{g-')) 

' V ' 

^loc 

= hoc{^,c{g-')) + lnp{ag),c{g-')). 

Note that / can be expressed in terms of (^, a) G x Og as 

/(^, a) = K{^) - Eiocia) - E^piCiga), a) = kod^, a) + /„p(C(fi'a), a), (6.13) 

where (7a G G is such that c{g~^) = a. This ga is determined only up to left multiplica- 
tion by Gg. Since Enp is Gg-invariant, the function a Enp{ga) is well-defined. Note 
that the Lagrangian of the strand (see (2.4), (2.13), and ( |2.14[ )) is exactly of the form 



(6.13), with ( = (^^{s, s'), K,{s, s')) G SE{3). Since a lnp{C{ga),o,) is a well-defined 
function of a G Cg one can ask why we insist in denoting /„p = lnp{C{ga),ci) instead 
of simply Z„p = Inpia) which is mathematically correct. The reason is that for the 



molecular strand we do not have an explicit expression for Inp : Og ^ M; see (2.14). 



detail in i.QA 



Note that (2.14) is exactly of the form l^p = lnp{C{ga) , ci) ■ This will be explained in 



Ellis et al. 



Dynamics of charged molecular strands 



46 



6.3.2 Recovering the modified Euler-Poincare approach 



By Theorem 6.3 we have seen that the Euler-Lagrange equations of a Gg-invariant 
Lagrangian Lq : TG M are equivalent to the affine Euler-Poincare equations for 
/ : X C»g ^ M, that is, 

^ - - ^ ^ • (6.14) 



dt 5i 



Recall that to write these equations, we need to extend / to g x V^*. Nevertheless, as 
we have shown, this extension does not affect the solution of these equations. For the 
molecular strand, there is an additional complication coming from the fact that the 
Lagrangian 

Ki^ o) = koci^, a) + InpiCiga), a) (6.15) 
being a well defined function of (^, a) G g x Oq, is not explicitly written in terms of a. 
Therefore, when computing the affine Euler-Poincare equations in concrete examples, 
there is still a dependence on Qa in the final equation, although we know that this 
dependence can be replaced by a dependence in a uniquely, by the results above. 
Let us apply the variational principle (6.11) to Lagrangian in ( [6.15^ . Let g{t) be a 
given curve in G. Take a family of curves gs(t) satisfying go{t) = g{t) and denote 
rj{t) := g~Ht)5g{t). Then 5 j!' /(^(t), c{g{t)'^))dt = implies 



d SI 



loc 



dt 6^ 



ad] 



Shoe _|_ S{lloc + I 



np J 



Sa 



o a 



dc' 



6 {lloc + I 



np J 



6a 



(6.16) 



However, Inp is 



Note that this equation is the abstract generalization of equations (3.29) and (3.32). 
Recall from the abstract theory that Inp depends only on a G Oq. 
given as a function of {({g),c{g~^)). Let 

61 



"np 



6a 



Tot 



denote the functional derivative of Inp viewed as a function of a G Oq only. Since 



every curve in Oq through a 
have 



c{g ^) G Cq is of the form c{g^ ^), where go = g, we 



d 

de 



Lp{C{9e),c{g^ )) 



e=0 



61 



np 



6a 
61 



Tot 



d 

de 



e=0 



np 



6a 



o a — dc^ 



Tot 



where rj := g ^6g. On the other hand 
d 

de 



61 



np 



6a 



61 



np 



6a 



Tot 



e=0 



inp{a9e),c{g:')) = { ^-j^Mav)) - l^-^.m + Mii) 



,r]a + dc{ri) 



Tot 



(6.17) 
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Equations (6.17) and (6.17) prove the following identity 



5a 



o a — do 



T 



Tot 



np 



Sa 



Tot 



6lnp 1 T ( ^^np 

oa — do ' 



Sa 



5a 



+ 9' 



-1 5lnp rp > 

5C 



equations (6.14) since 



where a = c{g ^). Using this identity in (6.16) we obtain the affine Euler-Poincare 

5lnp 51 



5l_ 
5a 



5a 



+ 



loc 



Tot 



5a 



Thus, the affine Euler-Poincare process recovers the results of the modified Euler- 



Poincare approach described in 53.2 



6.3.3 Hamiltonian approach 



We now explore the Hamiltonian counterpart of the theory, that is, the case of a 
Gjj^-invariant Hamiltonian '■ T*G M, defined only for a fixed value Oq G V* . 
As before, we do not suppose that Ha^ is induced from a G-invariant Hamiltonian on 
T*G ^V*. In particular, we do not know the expression of Ha for other choices of.. 
In particular, we do not know the expression of Ha for other choices of a. Such an 
HaQ is usually induced by a hyperregular G^^-invariant Lagrangian La^. 

As on the Lagrangian side, the reduced Hamiltonian is only defined on the submani- 
fold 

0* X C 5* 



and so Theorem 6^ cannot be applied. However, as is shown in the next theorem, 
the fact that the reduced motion is Hamiltonian on an affine coadjoint orbit remains 
true for this more general case. 

We need to introduce the affine coadjoint orbit O'^^^y The left V^* -valued group 
one-cocycle c : G ^ V* induces a left group one-cocycle a : S ^ (fl©^)* by 

a{g, u) = {uo c{g) - dc^{u), c{g)). 

The affine coadjoint action of S on s* is hence given by 



{g,u){fi,a) := Ad(^^„)-i(/i, a) + 



u 



The connected components of the coadjoint orbits (^O^^ 



UJ 



are the symplectic 

leaves of s* endowed with the affine Lie-Poisson bracket (6.7). Denote by S"^^^^ the 
isotropy group of the affine coadjoint action. 



Theorem 6.5 Let H, 



T*G 



be a Gl^- 



invariant Hamiltonian, where is a 



fixed element in V* . By Ga^-invariance, we obtain the reduced Hamiltonian h on 



0* X o:^ c s* 



defined by h{n,9g{ao)) = Ha^ig V)- 
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(i) Let a{t) e '^g{t)^ ^ solution of Hamilton's equations associated to Ha^ with 
initial condition fiQ G T*G = g* . Then {fi(t),a(t)) := {g(t)~^a(t),9g(^t)-^{cLo)) ^ 
s* is the integral curve of the Hamiltonian vector field Xh on the affine coadjoint 

orbit ao)'^ ) initial condition {fio,ao). Conversely, given jiQ G T^G, 

the solution a(t) of the Hamiltonian system associated to Ha^ is reconstructed 
from the solution a,{t)) of Xh G X {pi^^ ao)) ^^^^ initial condition {fio, oq) 

by setting a(t) = g(t)fi(t), where g{t) is the unique solution of the differential 
equation g{t) = g{t)sj^fj with initial condition g{Q) = e. 

(ii) Extending h arbitrarily to s*, Hamilton's equations on (^^(^/ioao)'^ ) ^^'^ 
written as 

where /^(O) = /iq, a(0) = Gq. 

Remcirk 6.6 It important to observe that the given Hamiltonian h is not defined 

on the whole dual Lie algebra s*. Part ii of the theorem states that the equations of 
motion can be nevertheless computed from the usual formula of an afiinc Lie-Poisson 
vector field by arbitrarily extending h to s*. Note that dh/d/j, and 5h/6a are only 
defined when one thinks of h as being defined on s*. It will be shown in the theorem 
that the extension of h does not matter. This difficulty will appear concretely when 
dealing with the molecular strand. 

Proof. 

(i) The action ^ of 5" on T*S induces an action of V given by 

{ah, {u,a)) {ah,v + u,a). 
Since y is a closed subgroup of S, this action admits a momentum map given by 

Jv{(^g, {u, a)) = a. 

Since V is an Abclian group, the coadjoint isotropy group of Qq G V* is 14g = V and 
the first reduced space {T*S)ao ~ J\/^('^o)/^ is symplectically diffeomorphic to the 
canonical symplectic manifold {T*G, flcan)- The action ^ of S* on T*S restricts to an 
action of G^^^V on Jy (do)- Passing to quotient spaces, this action induces an 
action of G^^ on {T*S)ao, which is readily seen to be the cotangent hfted action of G^^ 
on T*G. We denote by J^q : {T*S)aQ (flai,)* the associated equivariant momentum 
map, where g^^^ is the Lie algebra of G^^. Reducing (T*S)ao the point fj,ao := Ai|0a„, 
we get the second reduced space {{T*S)ao)^^^ = Ko{l^ao)/{Gl^)^^^, with symplectic 
form denoted by {^ao)nao- 
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By the Reduction by Stages Theorem for nonequivariant momentum maps [26], the 
second reduced space is symplectically diffeomorphic to the reduced space 

ao)/ 5'('^,ao)5 ^{^^,ao)) 

obtained by reducing T*S by the whole group S at the point (/i, ao) ^ s*. By 
affine Lie-Poisson reduction, this space is symplectically diffeomorphic to the affine 
coadjoint orbit 

endowed with the affine orbit symplectic symplectic form. 

Note finally that by the symplectic reduction theorem, any solution of Hamilton's 
equations associated to on T*G reduces to and is reconstructed from a solution 
of Hamilton's equations for the reduced Hamiltonian h^^^ : 3ao(^^ao)/i^ao)^^aa ~^ 
for a given momentum value fiao G (flao)*- have seen, this reduced space is 

symplectically diffeomorphic to the affine coadjoint orbit C s*, where G g* is 

such that fi\gc = ^ao- Thus, we can think of /i^^^ as being defined on O'^^^^y Viewed 
this way, /i^^^ is simply the restriction of the function h constructed from by 

hifi,Og{ao)) = H{g'^iJ,,aQ). 
Note that h is defined on any affine coadjoint orbit (9^^ with fixed oq € V* since 

(ii) We begin by recalling a general fact from the theory of Poisson manifolds. Let 
!f G C^lP), where P is a Poisson manifold and its Hamiltonian vector field. 
If L is a symplectic leaf of P, then X^\l = where the right hand side is the 

Hamiltonian vector field on the symplectic manifold L. In our case P = s* and 

Remark 6.7 (The case oq = and the charged strand) The Lagrangian 

Loivg) = K{vg) - EUc{g-^)) - E^Md) . c{g'^)) 



discussed in Remark 6.4 is hyperregular, thus it induces the GQ-invariant Hamiltonian 

Ho{ag) = K{a,) + ^/oc(c(r')) + Kp(C(^7), c(r ')) 
whose reduced expression on g* x Oq reads 

h{ii,c{g'^)) = ^^M"" + EioMg-^)) + E^p{a9).c{g-^))- 
As on the Lagrangian side, for (/i, a) G g* x Oq (or (yU, a) G 07 qJ, we can write 



K^^^ «) = + Eioc{a) + Enp{({ga), a), 



where ga E G is any group element satisfying c{g^ ^] 



a. 
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Remark 6.8 (AfRne coadjoint orbits and Noether's theorem) As we have seen, 
the solution (/i, a) evolves on an affine coadjoint orbit, for any G^^-invariant Hamilto- 
nian Ha^. If Lag is the Lagrangian of a simple mechanical system with symmetry then, 
by Noether's theorem, the solution (^, a) is constrained to evolve on the submanifolds 



6.4 Application to the charged strand 

In this subsection we apply the affine Euler-Poincare and Lie-Poisson reduction the- 
orems to the molecular strand. In order to give a more transparent vision of the un- 
derlying geometric structures, we consider the n-dimensional generalization described 



in Subsection 2.4, that is, we replace the interval [0, L] be an arbitrary manifold V 
and we replace SE{3) by the semidirect product S = 0(§)E of a. Lie group O with a 
left representation space E. Given a manifold V, we define the group G := J-^{T>, S) 
and the dual vector space V* := Q^(V,s) Q) J^(T>, E). The elements of the group G 
are denoted by (A,r), where A : V O and r : V ^ E. The elements of V* are 
denoted by {fl,T,p), where Q G Q^(V,o), F G Q^(V,E), and p:V^E. The space 
V* can be seen as the dual of \^ = X{V,s*) © J^{T>, E*), where X{V,s) is the space 
of s- valued vector fields on V. 
Consider the representation of G on V* defined by 

(A,r)(fi,r,p) = (Ad(A,,)(^],^),Ap) (6.19) 

where the adjoint action is that of S, acting here on functions defined on V, and 
Ap denotes the left representation of O on E, acting on functions. The main object 
for this approach is the group one-cocycle c appearing already implicitly in the 
definition of the variables fl,T,p in ( |2.11[ ), (2.20), and explicitly in (2.21). Recall 
that it is given by 

c(A, r) := ((A, r)d(A, r)"\ -r) . 

Let's verify the cocycle identity for the first component (A, r)d(A, r)~^. To simplify 
notation, denote Xi '■= (Aj, r^) G J-'(T>, S),i E {1; 2}. We have 

XiX2d(xiX2)"^ = XiX2d(x2 ^Xr^) = XiX2d{x2^)xT^ + X1X2X2 ^d(xr^) 
= Ad^, (X2dx2') +Xid(xr')- 

Since the second coordinate of ((Ai, ri)(A2, r2)) is equal to r = ri + Air2, we find 

c((Ai,ri)(A2,r2)) = 

= (Ad(Ai,ri) ((A2, r2)d(A2, r2y^) + (Ai, ri)d((Ai, r^)-^), -r^ - A^r2) 
= (Ad(A„.,) ((A2,r2)d(A2,r2)-i) ,-Air2) + ((A^, ri)d((Ai, n)'^), -n) 
= (Ai,ri)c(A2,r2) + c(Ai,ri). 
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This shows that c verifies the cocycle property (6.2) relative to the representation 
( |6l9| . ♦ 

Note that the first component of c is the left version of the cocycle appearing in the 
theory of complex fiuids; see [17]. Using the expressions 

{u, w, f) o {n, r, p) = (ad;^^ u' + w'oTi + fop, -n,w') , 

dc{uj, 7) = {—du, —d'j, —7), and dc^{u, w, f) = (div(u), div{w) — f) , 



the affine Euler-Poincare equations (6.5) become 

61 



(5.-ad:)- + (div-ad;^;--- 



5l_ 

51 51 



51 ^ 51 
07+^or + -op. 



(6.20) 



and the advection equations are 

dtfl + ad(^ Q = duj , 
(dt + cu) r = (d + fi) 7 
dtp + ujp = 'y . 



(6.21) 



Remark 6.9 To write these equations, we have supposed that the dynamics is de- 
scribed by a Lagrangian / given explicitly in terms of the variables (cu, 7, fi, F, p). 
Equivalently, we have assumed that / is induced by an affine left-invariant Lagrangian 
L defined on TG x V^*. As we have seen in §2.2[ such a hypothesis is not verified 
when nonlocal terms are taken into account. In this case, the affine Euler-Poincare 
and affine Lie-Poisson reductions are not applicable and one needs to restrict to a par- 



ticular value of the parameter oq, by using Theorems 6.3 and 6.5 For convenience, 
we first present the simpler case where the nonlocal terms are ignored. We shall call 
this case elastic filament dynamics for simplicity. 

6.4.1 Elastic filament dynamics and Kirchhoff's theory 

Suppose that the dynamics of the strand is described by a Lagrangian I = l{uj, 7, Q, F, p) 
defined ouqxV*, where = J^iV,s) and V* = n^{V,s) ®J^{V, E). The Lagrangian 
/ is induced by a left invariant Lagrangian L defined on TG x V* , where G = J-^{T>, S). 
Note that there is no restriction in the way I depends on the variables. In particular 
the dependence can be nonlocal. However, it is supposed here that / depends explicitly 
on the variables (cj, 7, il, F, p). Recall that such an hypothesis is verified for the 



molecular strand (2.4). 



Lagrangian of Kirchhoff's theory (1.11) but is not verified for the Lagrangian of the 
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The affine Euler-Poincare reduction applies as follows. Fix the initial values {Qq, Fq, po) 
and define the Lagrangian 

L{no,ro,po)i^^r) := L(A, r, fio, Tq, Po)- 

Consider a curve (A, r) E G and define the quantities 

{Q, F, p) = (A, rr\no, 70, Po) + c((A, r)-') 

= (AdA-i Qo, A"'(ro + fior), A" Vo) + (A'MA, A-'dr, A" V). 

and 

u = A^^A, 7 = A~V. 
Note that when the initial values no,Fo,po are zero, the definition of the variables 



a;,7,n,F,p coincide with those given in (2.11) and (2.20). 



Then the curve (A, r) is a solution of the Euler-Lagrange equations associated to 
L{no,ro,ro) on TG if and only if {u, 7, Q, F, p) is a solution of the Euler-Poincare equa- 
tions (le^tol). 



Of course, when V is the interval [0, L] and S is the semidirect product oi O = 5*0(3) 
with E = R^, then we recover from (|6.20l) the dynamical equation of the charged 



strand (4.12), since 



ad — — X and o — > x . . . 
These equations are the convective representation of Kirchhoff's equations. From 



(6.21) we recover the advection relations derived in Subsection 2.3 



6.4.2 The charged strand: general case 



Recall from ^2.2 that the Lagrangian of the molecular strand has the expression 

/ = kocii^, 7, ^, r, p) + inp{^, K, r). 



where linr is a local function of the form 



hocic^, 7, ft, r, p) = Kicj, 7) - Elocift, r, p) 



(6.22) 



and Inp is of the form 



l^p{^,K,r) = jj U{i{s,s'),n{s,s'),ns),T{s'))dsds', 

where 

U : SEiZ) xm? xm? and (^(s, s'), s')) := (A, r)~i(s)(A, r)(s'). 

Remark 6.10 (Two crucial observations) 
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1. The nonlocal Lagrangian l^p is induced by a S'0(3)-invariant potential Enp = 
Enp{A,r). Thus the total Lagrangian / can be seen as being induced by the 
5'0(3)-invariant Lagrangian Lq = Lo{A, A,r,r) given by 

Lo(A, A, r, r) = K{A, A, r, r) - Ei^c (c ((A, r)-')) - Kp(A, r), 

where K is the J-'(T>, S'i?(3))-left invariant extension of the kinetic energy K in 
(6.22). Note that we have replaced the dependence of Eioc on (f2,r,p) by a 
dependence on (A, r) through the cocycle c. The affine Euler-Poincare dynamics 
will yield the relation (f2,r,p) = c((A,r)^^) which allows us to recover the 
dependence of the potential on (f2, F, p). 

2. The group S0{3) is precisely the isotropy group 

= ^(P, 5E(3))^ = {(A, r) G G I c(A, r) = 0} 

of the affine action at zero. 

These two remarks allow us to obtain the dynamics of the molecular strand by the 
affine reduction processes described in Theorems |6.3 and 6^ As before, we choose 
to work with the general framework involving V and 0(§)E. The present approach 
is applicable to any (9-invariant Lagrangian 

Lo = Lo(A, A, r, r) : T[J^{V,0(i) E)] ^ R. 

Note there are no conditions on the dependence of Lq on the variables (A,r). In 
particular, Lq can be nonlocal, and may depend on the derivatives of A and r. An 
important class of such Lagrangians is given by 

Lo(A, A, r, r) = K{A, A, r, r) - P(A, r), 

where K is the kinetic energy associated to an (9-invariant metric on J-'{V,0(§) E) 
and the potential P is an O-invariant function on J^{V, 0(§)E). In particular, P can 
be nonlocal, or depend on derivatives of A and r; see (1.11) for an example. In the 
case of the molecular strand, K is assumed to be left-invariant and P is given by 

P(A,r) = Eio, (c ((A,r)-i)) + E„p(A,r), 

where 

E„p(A,r):= jj U {^{s,s'),K{s,s'),A~^dr{s),A^^dr{s'))dsds' 
ias,s'),Kis,s')) := (A,r)-i(s)(A,r)(s') G 0(S)E 
and one readily sees that Enp is (9-invariant. Recall that the cocycle is 

c((A,r)-i) = (A-MA,A-Mr,A-V) . 
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Thus, a straightforward and maybe useful generahzation of Enp is 

E„p(A, r) := [[ uUis, s'), k{s, s'), c ((A, r)-') (s), c ((A, r)-') (s')) dsds' . 



V 



Using Theorem 6.3 with Lq we obtain the same affine Euler-Poincare equations (6.20), 



where all derivatives are total derivatives. One can equivalently use the modified 
Euler-Poincare approach and obtain the equations 

(a,-ad:)|^ + (div-adj)^ 



51 61 ^ 61 

6^'''^6f'''^+6p''P' 



+ 



, dU , dU , dU , 

^">~Q^^^'^ ^) ~ - k{s, s') O — (s, s') 



/r. ^ 61 , r^. 61 61 

(9,+.)^ + (d,v = - 



+ 



dU 
dn 



s, s 



ds' 



ds . 



(6.23) 

Note that here the derivatives are not total derivatives, see the discussion in §6.3.2[ 
One can treat the Hamiltonian side in a similar way. As we have seen, the motion is 
Hamiltonian on affine coadjoint orbits. 



6.4.3 Conservation laws and spatial formulation 

In this paragraph, we generalize the approach of Section |4] and reformulate the equa- 
tions (6.20) for the generalized charged strand as a conservation law. We first need a 



n-dimensional generalization of formula (4.3). Given a Lie group G, a map g : V ^ G 



defined on a n-dimensional manifold V, s G V, and a g*-valued vector field w on V, 
we have 

Ad* [div(Ad*-iw)] =divM;-ad;^w;^=:div"w, a := g-^dg e n\V, g). (6.24) 



Using this formula, (4.3), the expression of ad* associated to the semidirect product 



0@E, and the equalities 

(a;,7) = (A,r)-i(A,f), (^], F) = (A, r)-M(A, r). 



we find 



Ad 



d_ 



61 61 

d f 61 61 

+ 



Ad 



dt \6uj^ 6^ 



ad* — + o — uj — 
^ 6u! ^ 57' ^7 



and 



Ad*^ div 



■. f 61 6l\ 

^. f6l 6l\ [ 61 ^ 61 61 
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Thus, equations (6.23) can be rewritten in the form of a conservation law, namely 

+ div Ad 



d_ 

Ft 



51 61 



(A,r)- 



51 SI 



Ad 



(A,r)- 



61 



51 



— op, — 
5p 6p 



(6.25) 



Using (6.1), the right hand side simplifies to 
61 



A A ^ /^A 

A— oAp +r o A — 
6p J \ 6p 



A 



61 
Tp 



A 



6l_ 
Tp 



0,A 



61 
Tp 



since p = A^^r. Note that this is the exact equivalent of the simplification (1.23) 
derived at the beginning of the paper. 

Such a conservation law is valid for each solution of the affine Euler-Poincare equation 
(6.1) associated to a GQ-invariant Lagrangian Lq : TG M. In particular, it is valid 



for the Kirchhoff's theory as we saw at end of 1.3.1 



A short computation shows that, in general, the previous conservation law reads 



d_ 
Ft 



+ dc^ g 



6l_ 
6a 



When qq is not necessarily zero, the previous formula becomes 



Ft 



Ad,-. ^ 



+ dc' g 



61 
' 6a 



Ad;-.(£o,-ao 



(6.26) 



(6.27) 



6.4.4 The fixed filament and its conservation law 



The equations (1.28) for a fixed filament can also be obtained by affine Euler-Poincare 



reduction. It suffices to apply Theorem 6.3 with the group G = J-'{V, O) 3 A, acting 
on the vector space o) x JF(D, E) 3 {Q, p) by the affine action 

{n,p) ^ eA{n,p) ■= (AdAi] + AdA'\Ap). 

Note that the cocycle is c(A) = (AdA^"'^,0). Using the expressions 

(M,/)o(fi,p) = ad;^^ u' + fop, 

dc{uj) = (— dco',0) and dc"^(M, /) = div(M), 
the affine Euler-Poincare equations ( 6.5[ ) become 

61 



op 



(6.28) 
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and the advection equations are 



dt^l + ad^ Q = du . 
dtp + up = . 



(6.29) 



Recall from ^ 1.3.2 that the Lagrangian for a fixed filament is of the form 

/ = lioc{cL>,fl) +lnp{^,p), 

hUu;, ft) = K{u:) -^j f{n{s))ds, P) = ' jj U{p{s),^{s, s'))dsds' 

where 

/ : ^ M, U -.R^x 50(3) ^ M, ^(s, s') := A-\s)A{s'). 

Using the relations uj = A^^A, Q = A^-'^A', and p = A~^Pq, where Po{s) := r{s) = 
(s,0,0)^, we obtain that / is induced by a 5'0(2)-invariant Lagrangian i^{o,r) = 
L(o,r)(A, A). Note that 5*0(2) is precisely the isotropy group of (0, r) relative to 
the affine action. 

These observations allow us to obtain the equations for the fixed filament by the affine 



reduction processes described in Theorems 6^ and 6^ Using the general framework 
involving V and 0(§)E, we obtain the equations 



(9.-adJ- + (div-ad,)- = -op + 



dU , dU I 
i{s, s')^{s\ s) - — (s, s)i{s\ s) 



ds' 



which coincides with (1.27) in the case of the fixed filament. Using total derivatives, 



these equations can be rewritten as (6.28) 



The general formula (6.27) yields the conservation law 
d 



dt 



Ad 



51 
5uj 



+ div 



Adt-i 



51 
5Q 



51 



AdI_M^op 



From the general theory it follows that the solution of the advection equations (6.29) 
in terms of A are given by f2 = A~^dA and p = A~^po- 

For the fixed filament, we choose V = [0,L], E = R^, O = 50(3), Pq{s) = r{s) = 
(s, 0, 0)^ and we get 



d_ 
dt 



51 
5u: 



+ 



d_ 
ds 



51 
In 



(6.30) 



Note that in this case, the torque does not vanish. The explanation is that the initial 



value Pq of p is not zero, so we need to use (6.27) instead of (6.26). .. Observe that 
we can write 

s 



Adl-i 1- X P = A— X Ap = A— X 
\5p J 5p I 
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More generally, the right hand side is 



A 



61 
'6p 



X r. 



where r describes the fixed filament. 



Note that the conservation law (6.30) does not appear in ^1.3.2 It is a particular 



case of the general formula (6.27). We believe that the derivation of this law through 



the affine Euler-Poincare theory is interesting and shows the breadth of application 
of our theories. 



7 New variables: Coordinate change and horizontal- 
vertical split 

In this section, we show that a drastic simplification of the equations arises under a 
particular change of variables. We shall assume that the Lagrangian / is only local. 
As far as we know, there is no general theory that can deal with the nonlocal term 
in the context of field theory, of which the present section is a forerunner; the field 
theoretic approach is developed for local Lagrangians in the next section. 
We first consider the case of strands. This change of variables will then be extended 
to the general setting of the previous section where [0, L] is replaced by a manifold V 
and SE(3) by an arbitrary semidirect product associated to a representation. 

7.1 Motivation in terms of covariant derivatives 



We can see from (1.1) that p, T, and 7 satisfy the following relations 

{ds + ftx)p = T, {dt + ujx)p = -f. 



Thus the reduced variables (1.1) lead naturally to two differential operators which 



can be interpreted as covariant derivatives, 

^^ = {ds + nx), -g = (9, + a;x). (7.1) 

With this interpretation we regard F and 7 as covariant tangent vectors above p, 

^ = T ^ = 7 (72) 
Ds ' Dt ' 



The operators from (7.1) also appear in the equations of motion (4.12) since we can 



write the second Euler-Poincare equation in the form 

D 51 D 51 51 



Dt5j ' Ds5r 5p °' ^''■^^ 
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When take (7.2) and (7.3) together we see that (7.3) is in the form of the Euler 



Lagrange equations where the partial derivatives have been replaced by covariant 
derivatives. With this interpretation in mind we can ask whether, by a change of 



variables, we can transform (7.3) to the canonical Euler-Lagrange form. In this section 



we find that such a change of variables does exist, and we give it explicitly. This line of 
enquiry leads us to consider in the subsequent Sections how the two sets of coordinates 
are related from a geometric point of view. 



7.2 The case of charged strands 

Consider the coordinate change 

J^([0,L],5o(3)) X J^([0,L],M3) X n\[0,L],5o{3)) x ^]l([0, L], M^) x J^([0, L], M=^) 

9 (lj, 7, n, r, p) ^ (p, p„ p„ Lj, n) e (7.4) 

J^([0,L],M3) X n^([0,L],M^) X J^([0,L],M3) x J^([0, L], so(3)) x n\[0, L],5o{3)) 



where we have defined two new variables 

p^ = T -n X p, p^ = j-uj X p. 



(7.5) 



We shall show that the equations of motion (3.3) and (3.4) have simple expressions if 



one uses horizontal and vertical coordinates. As far as we know, this transformation 
has not been noticed before, in either nonlocal or local setting. 

Notation 7.1 Assume that the Lagrangian / is local in the variables (w, 7, fl, F, p), 
that is, / = lioc- We shall denote by / the integrand of the Lagrangian / in terms of 



the new variables given by (7.4), that is, we have 



i{u;,j,n,r,p) 



7.3 Change of coordinates 

The action principle for the Lagrangian / = li^c yields 



= 5 / /(a;,7,f7,r,p)dt 



61 



61 



6Q, 



dt 



j KP,Ps,Pt,^,^)(isdt 



(7.6) 



61 



61 



61 



61 



61 
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Define free variations ^'(s) = A(s)^^5r(s) and = A(s)^^5A(s). As usual, \l/ 
denotes tlie antisymmetric matrix tliat is obtained from ^ by tlie hat map. Then, 
the following theorem holds. 

Theorem 7.1 The variations in Spg and 5pf- yield dynamical equations in the fol- 
lowing form 



61 SI 

{ds + nx) — + {dt + LJX) — 
oil 0(jJ 

SI ^ SI ^ SI 
Sp Spt Spg 







0. 



(7.7) 
(7.8) 



Remark 7.2 The derivatives in the equations (7.7) and (7.8) have now formally de- 



coupled, although the equations themselves must be solved simultaneously because 



the Lagrangian I depends on all the variables. Also note that equation (7.8) is equiva- 



lent, for local Lagrangians, to (3.4) with the covariant derivatives replaced by partial 



derivatives (but relative to the new variables). This gives a new interpretation to 



the right-hand side of (3.3) as being terms that arise from the induced covariant 
derivative. 



Proof. First, variations Sp^ and Sp^ are computed from (7.5) as follows: 

Spf = S^ — Suj X p — uj X Sp . 
ST — SVt X p — fl X Sp . 



(7.9) 
(7.10) 



Then, using the identities 



^7 
ST ■- 



Su: = Ti + u! X Ti , 
= S' + ri X S , 
Sp = -S X p + * , 
:* + a;x*-Sx7, 

*' + Ox*-sxr 



we find, for example, from the term involving the derivatives with respect to Pf , 



SI 



SI 



* + a;x*-Sx7 



~ {il + u xT^ X p — u X ^ 

d SI 



Exp 





SI 

-^^dtV^'Ji, 



7 X 



SI 
Spt 



P X w ) X — ,S 



(7.11) 
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where we have used the Jacobi identity simphfying two triple cross products. We now 
employ the .. We now employ the kinematic condition for the derivative of p, 

dtp = j-uj X p, 



to simplify the S term in (7.11) and obtain the following simple condition 



Analogously, 



61 



61 



Sp. 



d 61 
dt 6pt 

d 61 
ds 6p^ 



^ d 61 ^ 
^ dt6pt' 



d 61 ^ 
p X —- — , S 



ds 6p^ ' 



(7.12) 



(7.13) 



On completing the variational principle (7.6) for all variables, one sees that the only 
terms containing ^ are the derivatives with respect to p, p^ and p^. Due to (7.12) 



and (7.13), these remaining terms yield (7.8) 



On collecting the terms proportional to S, we notice another cancellation. As is 



evident already from (7.12) and (7.13), all the terms involving cross products with 



respect to p will cancel, as they will each be multiplied by the left hand side of (7.8) 
which vanishes. Thus, derivatives with respect to p, p^ and p^ will not contribute to 



the terms proportional to S, so that collecting those terms will yield exactly (7.7). 



There is another approach to performing the change of variables that highlights the 
decoupling. We key point is that we recognize two pieces of information we know 
about the variations 6p, 6pg and p^. First we consider the expression for 6p in terms 
of the free variations ^ and S. The relation is given by 

5p = * - S X p. 

This relation can be interpreted as saying that we can select any two of the variations 
S, and 6p as a free variation and the third variation is then determined. We find 
in practice that there are quantities such as 6fl that only depend on S. Therefore 
any selection of free variations must include S. This leaves us with a choice of ^ or 
6p as the choice for the second free variation. It is interesting to consider the choice 
of 6p. Indeed, since we have the relations 

Ps = dsP , Pt = dtp . 

we can express the variations 5pg and Pj in terms of our free variation 6p. 

6p^ = 6dsp = ds6p 

Similarly, 6pt = dt6p. Since (5f2 and 6uj only depend on S we have a complete 
description of the variations in terms of S and 6p which are given by 



6uj 



^Ps 



X S , 

= ds6p , 



6(jj - 
Spt 



s' + ri X s , 

^ dt6p, 



(7.14) 
(7.15) 
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which is obviously augmented by the trivial relation 6p = 6p. An alternative proof 



of Theorem 7. 1 can be given as follows 



Proof. Using the variations (7.14), (7.15) we obtain, for example, the following 



calculation in the variational principle. 

The terms arising from 6fl and 6u} are identical to before and only depend on I]. 
Therefore, we obtain the following equation from stationarity under the Xl variation, 

OUJ oil 

The second equation comes from terms proportional to 6p which is 

. SI SI SI 
Spt Sp, dp 



These are the required equations in Theorem 7.1 



Remark 7.2 Notice that this alternative proof does not require any cancelation of 
terms after the equations are derived. Thus, the variations do all the work for us. 
This opens up an interesting question. In some sense the choice of 6p as a free 
variation is optimal since no extra terms appear in the resulting equations of motion. 
We might also refer to the heavy top at this point and ask whether a similar change 
of variables might simplify the heavy top equations. The is negative 

but is nevertheless instructive. The crucial property that we used was to regard Sp 
as a free variation. Now, suppose we have an advected quantity a = A~^ao. This 
case appears in the heavy top as well as often occurring in fluid dynamics. Could we 
consider 6sl as a free variation? Unfortunately the variation 6sl is given by 

5a = — S X a . 

Therefore 5a is determined by S and we cannot interpret 6sl as a free variation. We 
shall investigate the geometric structure required for this approach in [8] 

Remark 7.3 This change of variables is not available in the classical Kirchhoff ap- 
proach because the variable p is absent in the classical approach. 

7.4 The general case 

We now generalize the previous results to the general situation described in Subsection 



2.4 Recall that in this case we have (A,r) G J^(V,S), {fl,T) G Q^(V,s), and p G 
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J^(V, E), where S = 0(§)E is the semidirect product of a Lie group O with a vector 
space E. 



Consider the variable p. Recall from (6.21) that we have the kinematic equation 

p = 'y — up . 

Assuming that the initial value of p is zero, we have 

dp = d(A- V) = -A^MAA- V + A-^dr = T-np. 

This motivates us to define the new variables ps G E) and pt G J-'{'D, E) which 

will play the role of space and time derivatives of p. They are naturally defined by 

Ps = T — Qp , and pt = '^ — ujp. (7-16) 

This change of variables defines a diffeomorphism from the variables (cj, 7, fi, F, p) to 



the variables {iu,il, ps, pt, p), and generalizes (7.5). In terms of the new variables, the 
local Lagrangian is denoted by / and we have 

i{p, ps,pt, t^, ^^)ds = /(cu, 7, n, r, p) . 



V 

There are two equivalent points of view to obtain the equations of motion in terms 
oil. 

The first one is to use a variational principle, as done before in the particular case 
of the charged strand. Using the constrained variations of to',7,f2,r,p given by the 
afiine Euler-Poincare principle, we obtain the constrained variations 

5uj = t + [cj, S] , 5VL = dT. + [Vl, S] , 

5pt = ^ - tip- T.pt , 5ps = d$ - dSp - Sp^ , 

and 

5p = $ - Sp. 

The second point of view is to compute the functional derivatives of / in terms of 
those of T. We find 

51 
ops 



and 

51 51 



51 
5uj 


51 

5u} ^ 


51 
5pt 


51 
M 


51 

~5n~ 




51 


51 


51 


51 




57 


5pt 


5f 


5ps ' 




51 


51 „ 


51 


51 



Ellis et al. 



Dynamics of charged molecular strands 



63 



These two ways lead to the same equations 

d 



dt 
d 51 



-ad* 



+ div^ — - 
5uj 5VL 



^ div — - — 

dt 5pt 5ps 5p 



(7.17) 



0. 



where div^ : X{V, o* 



J^{V,o*) is defined by 
div^ w := dww - ad^^ w' G J^{V, o* 



These equations coincide with (7.7) and (7.8) in the particular case V = [0,L] and 
S = SE{3). The other equations for the advected variables are computed to be 

Ps + ujps = dpt + iodp , 
Cl + adt^ Q = duj , 
P = Pt- 

We also know that dp = p^. Therefore, using the third equation, we obtain that the 
first equation is verified. Thus the last system can be replaced by 




8 The bundle covariant Lagrange- Poincare approach 

In this section we explain how the decoupled equations discussed above are covariant 
Lagrange-Poincare equations. The coordinate change will be interpreted as a trans- 
formation from the affine and modified Euler-Poincare perspectives to the covariant 
Lagrange-Poincare perspective. The corresponding Lagrange-Poincare equations are 
derived on principal fiber bundles by introducing a principal connection and splitting 
the configuration space into horizontal and vertical parts. Two equations occur, one 
horizontal and one vertical. In the case of strands we also have to take into account 
the continuous dependence of the variables on s. This leads us to consider a covariant 
version of the equations. We shall give various geometric structures that combine to 
give the required space. These geometric structures are introduced very effectively in 
the literature and the reviews of the various geometric objects are particularly based 
on [231 EH EE]. 

These field theoretical considerations work only for local Lagrangians. Since the 
classical infinite dimensional approach applies also to Lagrangians having a nonlocal 
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part, it is clear that an extension of the theory presented below exists also in the 
field theoretic framework. We defer to future work a development of the Lagrange- 
Poincare theory of Lagrangians depending on non-local variables. The local part of 
the equations of motion in the Lagrange-Poincare framework will, of course, have their 
non-local counterparts, equivalent to those derived in the Euler-Poincare framework. 



8.1 Covariant state space 

In this paragraph we shall introduce the covariant state space. The aim is to 
incorporate all the dynamical information into a single geometric object. 



We begin by noting that equations (4.12) have an exchange symmetry in their s and t 
dependences. Therefore, guided by the equations derived so far, we may treat s and t 
on an equal basis by introducing a spacetime, X := / x M. The dynamical quantities 
are then regarded as special vector bundle maps A : TX TSE{3). Such objects 
may be studied by considering the trivial fiber bundle 

Tixp : P := X X SE{3) ^ X, txxp{x, A,r) := x. 

The analogue of the state space TQ in field theory is the first jet bundle, J^P, of P. 

Definition 8.1 Given a locally trivial fiber bundle tcxp : P X , the first jet 

bundle tt : J^P ^ P of P is the affine bundle over P whose fiber at p E P is 

{j'P)p = {XeL (T,X, TpP) I TnxP o A = idr^x} , 
where L {T^X, TpP) denotes linear maps T^X TpP and x = ttxp{p)- 
The space J^P is called the covariant state space. 

Remark 8.1 It might, at first sight, seem unnatural to form the dynamics on a 
space of linear maps. After all, in canonical Lagrangian dynamics we consider tangent 
vectors, (g, q). However, in the canonical setting we could consider maps of the form 
Tq : TM TQ where the state space is TQ and M is time. Since Tq is a linear 
map on each fiber of TM we consider a basis of the image of each fiber given by 
T(^t,i)q ='■ {(l,<i), then T(^t,a)q = for all a G M, which is just rescaling of time 

viewed in a geometric way. When there is more than one independent variable we 
wish to capture the entire dynamics, independently of which direction is chosen in 
spacetime, therefore the notion of a linear map is the idea that generalizes most 
elegantly. In what follows one should think of jets as giving the 'velocities' in an 
arbitrary direction on spacetime. Thus it turns out that the first jet bundle, J^P, is 
a very natural state space since it is the analogue of the tangent bundle in the case 
where many independent variables are considered. 

In field theory one only uses certain sections of J^P, namely the holonomic or first 
jet extensions of sections of P. 
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Definition 8.2 Let a : X ^ P be a section of P, that is, vr^p o a = idx- The 
first jet extension of a is the map j^a : X J^P defined by j^cr(x) = T^a for all 

X e X. 

We see that j^a^x) G (J^P)o-(^), we differentiate the relation nxp o a = idx to find 

Ttxxp o Ta = id-Tx- This verifies that Ta G J^P. 

Given X = I x M. and P = X x SE{3), any section cr reads 

a{x) = {x,A{x),r{x)) G {x} x (50(3) ©M^), 

where x := {s,t). In this case we also have {J^P)a-{x) — L(TxX,T(^^x),r{x))SE{3)). 
Using this identification, we can write 

j^(j{x) = Tx {x, A, r) = (x, A(x), r(x), idr^x, A'{x)ds + A{x)dt, r'{x)ds + r{x)dt] . 



From (2.11) we conclude that the dependent variables that occur in the unreduced 
Euler-Lagrange dynamics are simply components of a first jet extension of a section 

of TTxp. 



8.2 Principal bundle structures 

Consider the natural principal 5'0(3)-bundle structure on SE{3) given by the pro- 
jection 

TisEis) ■■ SEi3) ^ W", nsEi3) (A, r) = A" V = p . 
The action of S0{3) is given by 

giA,r) = {gA,gr). (8.1) 

This principal bundle structure induces a principal S'0(3)-bundle structure on the 
trivial fiber bundle P = X x SE{3). We enforce the relationship ttxp {g ■ p) = '^xp (p) 
for all p E P. This relationship means that the group does not act on spacetime and 
this is reasonable because we do not want 5*0(3) to act on s or t in the applications. 
Explicitly, we have the following definition of the action of 50(3) on P: 

g{x,A,r) = {x,gA,gr) . 

Using this action we easily see that the reduced space is given by 

S := P/50(3) = (X X SE{3)) /50(3) = X x {SE{3)/SO{3)) =X xR\ 

The principal 50 (3)-structure on P is given by the projection induced on P by nsE{3)- 
That is, 

TTsp : P ^ S, vtep (x, a, r) = {x, nsEi3) (A, r)) = (x, p) . 
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Define the trivial fiber bundle ttxs : S — > X by 

TTXS {X, p) = X. 

To summarize, we are given a principal 5'0(3)-bundle structure TrsE{3) on SE{3) and 
a fiber bundle structure Hxp on P. From these we construct a new principal S0{3)- 
bundle structure, vr^p, on P and a new fiber bundle structure, ttxe, on E = P/S0{3). 
We note that .. We note that 



We summarize the considerations above in the diagram 



(8.2) 



S0{3) 



SE{3) 



SO {3) 



TTXP 



X 



id 



•""XE 



The section a{x) = {x, A{x), r{x)) of nxp : P ^ X induces the section x i— > {x, p{x)) 
of TTxE : E — > X. 

The tangent lift of the SO{3) action yields a free action on the jet bundle, J^P. 



9 



(x, A, r , id T^x.J^ds + kdt, r'ds + rdt^ 
= (x, gA, gr, id t,x, gA'ds + gAdt, gr'ds + grdt^ . 



The action is free because the action on 5*0(3) is free. Therefore we find that J^P 
is also a principal S'0(3)-bundle. In particular, J^P/ S0{3) is a manifold. We note 
that sometimes, for brevity, we omit the reference to x and idr^x in the explicit 
representation. 

Remark 8.2 At this point we could reduce by the 5*0(3) action on SE{3) to derive 
Euler-Poincare equations. The reduced variables are given as 

A-^r(A,r) = (e,A-V,A-^A'ds + A-^Adi,A-Vds + A-Vdt) 
= (p, Q.ds + ujdt, Yds + 'ydt) . 

This route is taken in the Euler-Poincare picture and results in the equations derived 
above. Again it is interesting to note that all of the dynamical quantities arise as 
components of a jet. 
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8.3 Principal Connection 

We introduce a principal connection that is needed to split J^P into horizontal and 
vertical parts and discuss its induced geometric structure. Recall that a principal 
connection on a principal G-bundle P is a g-valued one form on P that satisfies 

A{^p{p)) = t A{gv,) = Ad, A{p) , 

where gvp denotes the tangent lifted action of G on TP and is the infinitesimal 
generator 

d 



t=o 

For our particular 5'0(3)-bundle ttsp : P — > S, we make the choice 

A {x, A, r, v^, VA, u) = vaA~^ G 5o(3) , (8.3) 

for all va e TaS0{3), u eR^, e T^X. Note that this is the Maurer-Cartan con- 
nection for the structure group SO{?>). The choice of connection is actually arbitrary, 
but the particular choice above is well suited to the problem since it is not overly 
complicated. Also we recall that any vector Va G TaiS'0(3) can be written va = Arj 
where rj G 50 (3). The connection decomposes TP into the horizontal and vertical 
subbundles as follows: 

VerP = ker (Tttsp) = {{x, A,r]0, Ar], (AdAv) ■ r) \ rj e 5o(3)}, (8.4) 
HorA P = ker ^ = {{x, A, r; v^, 0,u)\ue TrR^, G T^X} . (8.5) 



8.4 Splitting TP/SO{S) 

In order to take advantage of the horizontal-vertical split of TP we give the induced 
global splitting of the vector bundle TP/S0{3) S. This is provided by the follow- 
ing vector bundle isomorphism, aA '■ TP/S0{3) — > TE ©s adP, given by: 




= Tn^pivp) ® [p, A(fp)]5o(3) ' '"p ^ TpP ' 



where adP := (P x so (3)) /SO {3) is the associated adjoint bundle to P. The 
quotient is taken relative to the left diagonal action given by 

(p,r/) ^ (hp,Adhr]), 

and the elements in the adjoint bundle are written \p,r]]so(3)- To check that aA is 
well defined we verify that 

TTrj:p{hvp) = T7rsp(fp), 
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and 

50(3) ■ 

To show that aA is an isomorphism we give its inverse, 



(nx,p) © [p,v]so(3)) = [^oT^inx,p)) + vpip)] 



SO{3) ' 



where p = {x,A,r) E P is such that 7rsp(p) = {x, p) and Hor^ denotes the hori- 
zontal lift of V(^x,p) = {x, p,Vx,u) G T(2. p)S to TpP with respect to A. It is given 
by 

(x, p, v^, u) = (x, A, r, t;^, 0, Aw) . 



Hor(^,Y,r)^ 



Remark 8.3 (The choice of connection) As we have seen in (8.3), a natural choice 



of connection is the Maurer-Cartan form for the structure group, dAA ^. 



8.5 Properties of adP 

We shall need various properties of the adjoint bundle to derive the Lagrange-Poincare 
equations; we review them here. 

We can give adP a Lie algebra structure on each fiber. The vector space structure is 
given by 

[P' V]sO{3) + « [P' ^]sO{3) = [P^V + ai^]sO{3) 

and the Lie bracket is given by 



[p,v] 



SO(3) 



[P, ^] 



SO{3) 



[P^ [^'H]sO(3) 



The principal connection A induces an affine connection on the adjoint bundle adP. 
It is known that the covariant derivative of this affine connection is given by 



£1 
Dt 



SO{3) SO{3) 



(8.6) 



(see, for example, [27j, Lemma 2.3.4). This formula allows us to define a covariant 
derivative of any section ( : X adP. Note that the adjoint bundle adP has base 
S. We view it now as a bundle over X. It is important to note that the composite 
bundle adP — ^ X is not a vector bundle, in general. The covariant derivative of the 
section ( is defined by using the formula for the covariant derivative of adP S 
induced by the principal connection A on vrgp : P ^ S. We define 



Dt 



(Coc)(r), 



(8.7) 



r=0 



where c(r) is a smooth curve in X such that c(0) = x and c(0) = U E T^X. Con- 



cretely, denoting ({x) = [p{x),T]{x)]so{3)i formula (8.6) gives 

V^C(x) = [pix),dvix)iU) - [AiTMU)),vix)]]so(3) ■ 
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Let us note that the vector bundle adP — > S is in our case trivial Indeed, the map 

[(x, A, r), r/]so(3) ^ {{x, A"V), AdA-i r]) (8.9) 
is a vector bundle isomorphism from adP to S x so (3). In this trivialization, using 



the connection (8.3), the formula for the covariant derivative (8.6) becomes 

^(x(r),p(r),e(r)) = (x(r),p(r),e(r)). 



Similarly, if ?7 G T^X, formula (8.8) becomes 



V^C{x) = {x,pix),d^{x)iU)), where = (x, p(x), ^(x)). (8.10) 

Had the bundles been nontrivial, the formulas for the covariant derivatives would be 
more involved. 

8.6 Splitting J^P/SO{3) 

Having introduced the connection that splits TP/ 5*0(3) we now wish to use it to split 
the reduced covariant state space J^P/ S0{3). This is easily achieved by regarding 
the jets as linear maps and composing with a^. Therefore we split J^P/ S0{3) by 
splitting the image of the jets in TP/ S0{3). 

Consider a section a of the fiber bundle nxp '■ P ^ X and its first jet extension 
jV(x) = T^a e (Jip)^(^). We compose [Ta]so{3) ■ TX TP/S0{3) with the 
vector bundle isomorphism aA '■ TP/S0{3) TS ©s adP over S and obtain the 
following equality in the fiber over vrsp(o"(a;)): 

aA o [Txcr]so(3) = {Ta{x)T^T.p o T^a) ® [cr(x), A o T^cr] (.^(g) 
= {t^t.p oa)®[a{x),Ao T^-cr] ^^^(s) • 

Using TTxE o TTsp = TTxp, we have 

TTxs o (ttep o a) = -nxp o a = idx--- 
This shows that ttsp o cr is a section of the fiber bundle ttxs : S ^ X. If we denote 

CTi = ttep o cr, A (vp) := [p, A (t^p)]5o(3) ' ^^{x) := A o T^a, 

for all Vp G TpP, then the reduced jet [j^<j{x)]so{3) ^ (J^P)/S'0(3) is expressed as 

aA o [jV(x)]5o{3) = "A o [Ta;Cr]<jo(3) = T^ai ®Ao T^a = T^p © cr2(x), 
since Ui = /?. Note that this element lies in the fiber 

{j'^U^x) X L(T.X,(adP),^(,) 
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over cTi(x) = ttxs{o-{x)) G S. In particular, there is a fiber bundle isomorphism 

J^P/S0{3) ^ J^S XsL(TX,adP) 

over S. Using the equality a{x) = (x, A(x), r(x)), the exphcit description of the 
quantities appearing in the reduced jet are: 



(^2 



(x,A-V(x),r,(A-V)), 

{x, p{x), p^{x)ds + p^{x)dt) e (J^S)<^(^), 

[{x,A{x),r{x)),T,AA-' 

((a;,A-V(x)),A-XA), 
{x, p{x),fl{x)ds + u){x)dt) , 



by the relations (2.11), where = denotes here the vector bundle isomorphism (8.9). 
Thus, the reduced jet [j^cr{x)]so(3) associated to a{x) = (x, A(x), r(x)) is represented 



in the trivialization (8.9) by 



Tx.p©a2(x) = (x,A-V(x),T,.(A-V),A-iT,A) 

= (x, p{x), Pg{x)ds + p^{x)dt, fl{x)ds + uj{x)dt) . 



(8.11] 



Therefore we have recovered the new coordinates given in (7.4) and the reduced 
Lagrangian reads 

i{Tp, (T2) = i{p, p„ pt, uj, n). 

We summarize the spaces involved in the previous discussion in the following diagram: 
adSE{3) adP 



TSE{3)/SO{3) 



TP/S0{3) '^--^/'''^^\ TX 



Ttt^p/SO{3) 

^ TE 



id 
TX. 



We recall that the projections ttxp and vrsp are S'0(3)-invariant, therefore they nat- 
urally induce projections tit.p / S0{3) : P/S0{3) — > S and t^xp / S0{3) : P/S0{3) — >• 
X, respectively. The relationship between the variables in the affine Euler-Poincare 
and covariant Lagrange-Poincare equations can be illustrated in the following dia- 
gram: 
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(p, nds + wrfi)so(3) ^ L{TX, adSE{3)) 



L{TX, adP) 



(p, Yds + jdt, nds + ujdt) e L{TX, TSE{3) /S0{3)) 

Tt:se(3}/SO{3) 

(p, p,ds + ptdt) e L{TX, TM3) 



J^P/SO{'S) L{TX,TX) 



Tttsp/SO{3) 

> jiE 



Tttxi 



L{TX,TX). 



The affine Euler-Poincare variables [j^o']soi3) — {p,^ds + ujdt,Tds + jdt) appear 
in the middle horizontal sequence whereas the covariant Lagrange-Poincare variables 
(p, Pgds + Ptdt, fids + ijjdt) appear in the top and bottom rows. 

8.7 Reduced Variations 

Let 0" : X ^ P be a section of the fiber bundle ttxp : P ^ X. If (Xg : X ^ P is a 
curve of sections with ctq = o", that is, a^^x) = {x, Ae{x),re{x)), Aq = A, and Vq = r, 
define the variation 



6a{x) 



d 
de 



{x,K{x),r,{x)) e T^^^)P. 



e=0 



Splitting Sa{x) into its vertical and horizontal parts relative to the connection A in 
the principal S'0(3)-bundle vr^p : P ^ S (see (8.4), ( 8.5[ )), we obtain 



6a{x) = (x. A, r, 0, 6A, 5r) 

= {x, A, r, 0, 5 A, 5AA- V) + {x, A, r, 0, 0, 5r - 5AA-^r) e T^(^)P. 

To compute the vertical variation of [j^o']soi3), we consider curves that perturb 
ao = a : X P along the group orbits, that is, 

<7e{x) := exp(£:^(x)) ■ a{x) = {x,exp{e^{x))A{x),exp{e^{x))r{x)) 



where ^ : X so(3). By (8.11), in the trivialization (8.9), 

[jVe(x)]^^(3) ^ (x,A-V(x);T,(A-V),(exp(£e(x))A(x))-'T,.(exp(£OA)) • 
Taking the 5-derivative of the right hand side we get 



aix 



d 



e=0 



SO(3) 



50(3) • 

= (^x, A-V(x); 0, (S' + ri X S) rfs + (^S + w X S j dt^ , (8.12) 
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where := AdA(x)-i ^{x). 

To compute the horizontal variation of [j^cr]soi3), we consider curves that perturb 
(To = 0" : X — i> P such that Sa is horizontal. In view of (8.5), a curve giving a 



horizontal 6a is a^ix) := (x, A(a;), rj(x)). Therefore, for such a curve we get 



50(3) 



(x, A-V,(x); r,(A-V,), A(x)"ir,A) 



and hence 

5^ [M 



50(3) 



e=0 



SO(3) 



(x,A-V(x);T,(A-i5r),0) (8.13) 



(x, p(x); T,(5p), 0) = (x, p(x); (5p)t + (M. ds, 0) (8.14) 



Remark 8.3 The free variations 6p and S are now recognized as being horizontal 
and vertical variations. This is the reason for the decoupled form of the resulting 
equations. If we had defects in the strand and therefore our connection had non-zero 
curvature, then the equations would not decouple completely. 

8.8 Variational Principle 

Having derived the reduced horizontal and vertical variations, we may now derive 
the horizontal and vertical Lagrange-Poincare equations. For the vertical variations. 



using (8.12) and 



we obtain 
d 



de 



X 



e=0 



Xr(bv.] 



51 
6a2 



50(3) 



p, —dt + —ds 



6uj 



6n 



30(3) 



dx 



P, —dt + ir^Os 



sn 



SO{3) 



X 



51 



p, ((9t + LJx) — + ((9s + rix, 

ou> oil 



p, ( S + X S ) rft + (S' + X S) 

51 



)dx 



SO{3) 



SO{3) 



SO{3) 



)dx = 0. 



Therefore, the vertical covariant Lagrange-Poincare equation is 

ou: 511 



Similarly, we derive the variational principal for horizontal variations and obtain. 
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using (|8.13l), 



d_ 

de 



=0 Jx ^ 



X 



X 



61 



5p) + 



61 ^ 61 
7 9t — 



50(3)) 


dx = 




61 


Sps 


Sps 


,5P^ 



ds , 6pf-dt + 6pgds )dx 



X 



Therefore, the horizontal covariant Lagrange-P oincare equation is 

1 ^t-r = 0. 

op 6pt 6p^ 

Upon putting these together, we find that the covariant Lagrange- Poincare equations 
are 



\ \ 

{dt + u:x) — + {ds + nx) — 

ouj 611 

61 ^ 61 ^ 61 

X 

op 6pt 6p, 








(8.15) 
(8.16) 



8.9 A circulation theorem 

The Kelvin- Noether Theorem tells us about the solutions to the Euler-Poincare equa- 
tions in continuum mechanics. There is an analogue in the covariant picture that is 
described in this section. Denoting by div^ the divergence relative to the variable 

X = (s, t) e [0, L] X M, we have 



61 



61 



6^ 



6u! 



div, Adt-i —ds + Ad:_i —dt ]=ds{ Adt -1 T7^]+dt{ Adt -1 



61 \ 



6n J 



61 
6uj 



61 



61 



61 



AdA-i 9s— + X — + dt— + X 



'6n 



6n 



6(jj 



61 
6u} 



by (8.15). Using the divergence theorem, we find 




div^ ( Ad^-i --ds + Ad^-i ^dt ] dsdt 
Oil duj 



as 



Adl-i ^ds + AdU ^dt]-ndi 
Oil Ou: 



as 



61 , 61 , 
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where n is the outward pointing unit normal to the boundary dS and we used the 
identity 

ds + ^dt) -ndi = ^dt - ^ds. (8.17) 



6fl 6u) J 6fl 6u) 

Thus, we obtain the circulation theorem 



(61 61 \ 

— dt- — rfs = 0. (8.18) 
6il 6u J 

8.10 Generalizations of the molecular strand 

Essentially the problem of the molecular strand is a complex filament. We have a 
filament in and also specify that the filament has microstructure that we describe 
by attaching a group element to each point of the filament. This setup can be general- 
ized in various interesting ways. First we can consider the multidimensional problem, 
that is, complex sheets and related structures. A molecular sheet is an object where 
molecules are bound in a two dimensional surface. In this case we can still con- 
sider RCCs of charges represented by 50(3), but this time the S'0(3) quantities are 
attached to points on an embedding of M? in M^. We could consider another gener- 
alization. Instead of having a filament with RCCs described by 5*0(3), some other 
group could describe the microstructure. For example, the group 50(2) was used in 
[13]. We could even extend the group to take account of quantum phenomena. We 
might want to make both of the above generalizations and consider, for example, the 
spin sheet. 

Motivated by these considerations, we quickly indicate here how to generalize the 



covariant Lagrange-Poincare approach to the setting of Subsection |2.4[ that is, the 

case of n-dimensional strand with an arbitrary Lie group structure O. 

Consider the [n + l)-dimensional spacetime X := P x M and the trivial fiber bundle 

■Kxp- P ■■= X X S ^ X , 

where S = 0(^E. A section a of P reads 

cr(x) = (x, A(x), r(x)), X = (s, t) G X , 

and its first jet extension is 

J V(x) = (T,A, T,r) = (dA(x) + k{x)dt, dr{x) + rix)dt) , 

where d is the partial derivative with respect to space (that is, the derivative on D), 
and the dot is the partial derivative with respect to time. 
There is a natural (9-principal bundle structure on 5 given by 

T^ES ■■ S E, TlEsi^, = A" V = p. 
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This principal bundle structure on the fiber S, induces a principal O-bundle structure 
on P given by 

TTsp : P ^ X X 7rEp(x, A, r) = (x, A" V). 
There is a natural connection A on ttsp : P ^ := X x E given by 

which allows us to identify the reduced jet bundle J^P/O with the fiber bundle 
J^S L{TX, adP). Using the same notations as before, we have 

"A o [j^cr{x)]o = T^p ®Ao T^a = (x, p(x), dp(x) + p{x),Q{x) + uj{x)dt), 



by (2.20). The vertical and horizontal variations being given by 

S" [f(r{x)]^ = (x,A-V(x);0,dS + + (S + [u,E])dt), 

d"" [jV(x)]^ = (x,p(x);d(5p) + (5p)irft,0), 



we find that the vertical and horizontal Lagrange-Poincare equations are 

0, 
0. 



61 



{dt - ad* ) ^ + (div - ad^; 



61 



61 
61 



6p 6p^ 6p, 



Of course, as expected, these equations coincide with equations (7.17) obtained from 
the affine Euler-Poincare equations (6.20) by the change of variables (7.16). In Sec- 



tion |9j we will obtain the same equations by applying directly the theory in [23], using 
the group structure of S. Note that the approach we have used here does not use 
the group structure of S, and is expected to be applicable to more general situations 
such as the molecular strand on the sphere. This is explained in the next subsection. 
We have 



div^ ( Ad^-i ^ + Ad^-i ^dt 



div ( Adl-i 



6r\ 

6q) 



+ dt AdU 



61 \ 



- 



6iu J 



Adt-i div 



61 



adr 



61 



61 



ad*, 



61 



6n ^ ^ 6n ' "'6uo ^""^ 6uj 
Using the divergence theorem, we obtain the zero fiux theorem 



0. 



dV 



A SI . _ 61 ^ 

Ad^-i + Ad^-i —dt 



6n 



6uj 



■nda = 0, 



where n is the outward pointing unit normal to the boundary dV of a given domain 
V (Z V X K. and da is the induced boundary volume element of dV. 
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Remark 8.4 We remark at this point that we have purposefully not referred to 
SE{3) in this section. Our entire method of reduction and derivation of (7.7), (7.8) 
makes no use of the group structure of SE{3). This is the first application we know 
that uses this method of reduction. The closest theory that one finds previously is in 
[23] where covariant Lagrange-Poincare equations are derived in the case in which a 
principal fiber bundle is reduced by a subgroup of the structure group. Naturally that 
theory is related to this approach in which we imposed less restrictive assumptions 
on P and renounced the SE{3) group structure. Of course, since in this case P is 
also a principal S'i?(3)-bundle, the charged strand is in the intersection of these two 
reduction theories. A treatment of the charged strand by a direct application of [23] is 
presented in the next section. The approach used here may be regarded as a covariant 
formulation of the Lagrange-Poincare equations derived in [27]. It also answers the 
call for the approach used in [23] to be applied to general field theories. One finds 
in general that under certain equivariance conditions on the trivializing maps of the 
bundle vr^p, a variational approach to field theories may be formulated that reduces to 
the case considered in [23] when the fibers of ttxp are Lie groups. This new reduction 
procedure also applies to a wider category of fiber bundles which includes the case 
where ttxp is trivial, but Q is any G-bundle, as often occurs in applications. This 
generalization can be seen in the present application by varying our assumptions. We 
note that the method above is capable of dealing with the problem of the molecular 
strand on a substrate or on any submanifold of M^, whereas the method developed in 
[23] does not treat that case. Of course the application of the above theory to such 
a system is subject to symmetry breaking. This is an interesting direction for future 
research, but it is beyond the scope of our present considerations. 



9 The subgroup covariant Lagrange-Poincare ap- 
proach 

In this section we allow ourselves the group structure of SE{3) and apply the results 
formulated in [23] to the charged strand. More precisely, we see the principal bundle 

nsEis) : SE{3) ^ M^ 7rsi<;(3)(A, r) = A'^r = p, 

as being associated to the left subgroup action of S0{3) = 5*0(3) x {0} on SE{3). 
Using the composition law in SE{3) 

(Al, ri)(A2, ra) = (A1A2, ri + Aira), 

we obtain that the subgroup action is given by 

S0{3) X 5^(3) ^ SE{3), Ai(A2,r2) = (AiA2,Air2). 



We thus have recovered the action (8.1 ), and the projection 7r5£;(3) identifies an equiv- 
alence class [A,r] in SE{3)/SO{3) with the vector A'V e M?. 
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Therefore, we can obtain the equation of the molecular strand by reducing the prin- 
cipal S'£'(3)-bundle P by the subgroup SO (3). Such a theory is developed in [23] . 
and is applied below directly to the n-dimensional generalization of the molecular 
strand. The difference with the approach described in the previous Section lies in the 
fact that Section [s] does not use the group structure of SE{3) and the fact that the 
principal bundle nsE{3) '■ SE{3) M.^ is associated to a subgroup action. 

Consider the manifold X = VxK. and the trivial principal S-bundle, P = X x S X , 
where 5* is a semidirect product group 0(^E. Since P is trivial, the first jet bundle 
is given by J^P(x,g) = L(TxX,TgS). A section a of P reads 

a{x) = {x, A{x), r{x)) , x = {s,t) , 

and its first jet extension is 

jV(a;) = {x,idTxX,Hx),r{x),T^A,T^r) 

X, id-xxx, A(x), r(a;), dA + Adt, dr + rdt 



where d denotes the derivative with respect to space and the dot denotes the time 
derivative. 

We can see P as a O-principal bundle over S := P/O = X x E, relative to the 
projection l : {x,A,r) t— > (x,A~^r). Suppose that we have a (9-invariant Lagrangian 
density C defined on J^P. This Lagrangian induces a reduced Lagrangian density 
/ : J^P/O M. On the principal bundle P P/O we consider the fiat principal 
connection 

where Vx G T^X, v\ G T\0 and G T^E. Using this connection, we have the fiber 
bundle isomorphism J^P/O = J^{P/0) x^ L(TX,adP) over S. Note that in our 
particular case, the vector bundle adP ^ S is trivial, see ( |8.9 ), and can be identified 



with S X 0. Moreover, the connection A is identified with the trivial connection, see 



ds.ioD . 

We recall now from [23] the covariant Lagrange-Poincare reduction, adapted here to 
the case of a semidirect product S = O E, and to the fact that P, as a 5'-principal 
bundle, is trivial. 

Given a section a = (A, r) of P ^ X, we consider section cti of P/O X defined by 
(7i{x) := L o cr(x) = A{x)^^r{x) = p{x), and the section o"2 of L{TX, o) — > X defined 
by (72 (x) = A~^Tj;A = i7 + udt. The following are equivalent. 

• 0" is a critical point for the variational principle 

6 [ C{j'a) = 0; 
J X 

• a satisfies the Euler-Lagrange equations for £; 
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the variational principle 

6 [ [(jVi,a2)=0 
Jx 

holds for arbitrary variations 6ai and variations of the form 

Sa2 = d^r7+ [a2,r]] , 
where r] is an arbitrary section of L{TX, o) — > X; 
the sections ai,a2 satisfy the covariant Lagrange- Poincare equations 



61 
6ai 



div^ 



div^^ 



61 \ 
61 



(9.1) 



^' 6ao """^^^ 6a. 



where div;^ denotes the covariant divergence associated to A and acting on 
X(X, 0*). Note that here div;^ = div^. 

Using the decomposition X = P x M, we can write 

/:(jV) =/:(A,dA,r,dr) 

and 



Hence, we obtain the equality 



61 
6a2 



61 61 
6VL 6uj 



a. 



Since div^ '^i^Tt ^"^^ ^"^^x '^2 = ^iv, the second equation of (9.1) reads 



d 61 61 61 61 



The first equation reads 



61 d 61 , 61 

div - — 

6p dt6pt 6ps 



0. 



We have thus obtained equations (7.17) by covariant Lagrange-Poincare reduction. 
Of course, when 0®E = 50(3) ©M^ = SE{3) and V = [0,L] we recover equations 
(8.15) and (8.16) for the molecular strand in the new variables. 
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10 Formulation of nonlocal exact geometric rods 
in terms of quaternions 

Quaternions allow for a simple, elegant and useful method of describing the local 
orientation of a curve. It is thus natural to seek a representation of our derivation in 
previous sections that expresses the strand equations in terms of quaternions. The 
quaternion representation is natural, for example, in formulating the equations of 
motion for elastic rods in terms of the corresponding Euler parameters. As far as we 
are aware, a treatment of continuum rod theory in terms of quaternions in the nonlocal 
sense presented here does not appear in the literature. We shall see how the nonlocal 



contribution (3.30) appears as an imaginary part of a certain quaternion, thereby 
making the connection to other work. This is accomplished by mapping quaternions 
(elements of SU{2)) that describe rotations into purely imaginary quaternions, or 
vectors, that are elements of su(2) ~ so (3) ~ M^. 

Remark 10.1 This section simplifies the formulas and avoids extra factors in the 
integrals by assuming the strand to be inextensible; so that |r(s)| is identically equal 



to unity and the parameter s is the arc length. See (3.29) and (3.32) for the case 

|r(.)|^i. 

Let us associate a quaternion q = (go, q) with every point on the curve s. That 
quaternion describes the local rotation of an orthogonal frame if the condition ||q{| = 
Q'o + I^P = 1 is satisfied. Then, go = cos a/2, with rotation angle a and q = sin(a/2)n, 
where v is a unit vector around which the axis is rotated. These are the Cayley-Klein 
parameters of the rotation 



Remark 10.2 To simplify the notation, 
quaternions, considering them as vectors. 

b = 



we use bold symbols for purely imaginary 
For example, if q is unit quaternion, then 

qaq* 



means 

(0,b) = q(0,a) q*. 

As before, we assume that the interaction potential depends on the distances between 
point charges that are attached to each point r{s,t) by rigid rods of the length 
r7j(s). The new position of the charges will be qr/^q*, where denotes the quaternion 
multiplication and q* is the quaternionis conjugate of q. The point charges are then 
positioned at the coordinates in real space r{s,t) + c\r]i(]*, and the distance between 
point charges is thei]]^ 

dkAs, s') = \ris) - ris') + q{s)rj,{s)q* (s) - q(s')r7^(s')q*(s')| • (10.1) 
■^AU these variables depend on time t as well as s, but the time variable t is suppressed. 
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This is simply (2.3) written now in its quaternionic form. Following (2.6), we perform 
the Lie-Poisson reduction as follows (remember that ||q|| = 1 and qq* = c, where 
c = (1,0) is the unit quaternion): 

4,m(s, s') = \r{s') - r{s) + q{s)r],^{s')q*{s') - q{s)r]^{s)q*{s)\ 

= \q*{s) (0, r{s') - r{s) + q{s')rj,{s')q*{s') - q{s)rj^{s)q*{s)) q{s)\ 

= l3(s, s')p{s')f{s, s') - p{s) + 3(s, s')r7fc(s')3*(s, s') - r}^{s)\ 

= \k{s,s') + ri,{s) - ^is,s')ri^{s')f{s,s')\ , (10.2) 

where ^(s, s') = q*(s)q(s') is the coupling between the frames and the quantity 

p{s) = q*(s)r(s)q(s), 

is the distance vector connecting the points r{s) and r{s') transformed according to 
the inverse rotation of the frame at the point s. We have also defined 

A€(s, s') = 3(s, s')p{s')f{s, s') - p{s) . (10.3) 

The nonlocal part of the reduced Lagrangian depends on the variables p and 3. The 
local part which describes elastic deformation and inertia, can be reduced to functions 
of m = q*q' and D = q*q where the prime denotes the derivative with respect to s and 
the dot is the derivative with respect to t. The quaternions m and belong to the 
Lie algebra of the Lie group of all unit quaternions, and is isomorphic to the space of 
purely imaginary quaternions, or vectors, as can be seen by differentiating q*q = e. 
The commutator is then mapped into twice the vector product of the imaginary parts 
of the quaternions. 

We can again split the reduced Lagrangian It in the local and nonlocal parts 
1 = y"//oc(o(s),m(s),7,r,p)ds + j j U {n{s,s'),i{s,s'))dsds' ■=lioc + lnp- (10.4) 



Here, 7, F, p are defined as in (2.11 ). The equations of motions then follow from the 



minimization of the reduced action 



5S = 5 j /(D(s),m(s),p(s),K(s,s'),3(s,s'))dsds'dt = 



;io.5) 



Here D and s are elements of Lie algebra of quaternions with fixed absolute value. 
They are purely imaginary quaternions, = (0,ci;/2) and s = (0,f2/2), which are 
isomorphic to vectors in M^. The factor of 1/2 is necessary for lo and f2 to be exactly 
the vector angular velocity and strain rate, respectively, in correspondence with the 



previous section. Thus, we can re-write (10.5) using vector quantities instead of 
quaternions whenever possible: 



6S = 6 Inp (cx;, 7, r, p) dsdt + 6 U (k(s, s'), 3(s, s')) dsds'dt = . (10.6) 
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We obtain 
SS = 

+ 



61 



loc 



6n 



6n ) + 



61 



loc 



6uj 



,6uj) + 



61 



loc 



6p 



61 



loc 



6j 



, 6j 



61 



loc 



6T 



6T^dsdt + 




61 



np 



Si 



61 



np 



6k, 



6k) dsds'dt = 0. (10.7) 

arly to 
(10.8) 

(10.9) 



If we now define s = (]*6q = (]^^6q as the free variation in q, we obtain, similarly to 
previous section: 

5d = OS — SO + 5 = [d , s] + S , 

for the time derivative and 

5m = ms — sm + s' = [m , s] + s' , 

for the space derivative. Note that since q*q = 1, 

q*q + q*q = q*ci + (q*q)* = 2Ret) = 

and, analogously. Rem = 0, which means that and m are purely imaginary quater- 
nions, or vectors. This gives the variation of the first two terms in (10.7). We now 
remember that Re = and Re s = since they are elements of the corresponding 
Lie algebra, so = (0,a;/2), m = (0, n/2) and s = (0, S/2). Then, ( flasj ) and ( [lO^ 
become vector equations: 



and 



(5o = (0, (5cx;) = (^0 , X S + S 
5m= (0,517) = (0, X S + S') . 



^10.10) 



'10.111 



Now, computations of the first three variations in (10.7) can be done analogously to 
those in the previous section, as they involve vector quantities. The only exception 
is the computation of 6p in the nonlocal term as it must be computed in terms of 
quaternions. We have 



6p = 6{q*{s)r{s)q) 

= -q*6qp{s) + p{s)q*6q + q*{s)6r{s)q{s) = 2p x S 

where we have defined the free variation 

*(s) := q*{s)6r{s)q{s) . 



10.12) 



(10.13) 



and used the fact that for purely imaginary s = (0, S), sp — ps = p x S. Then, the 
variation we need to compute the variation of the nonlocal part of k(s, s') as follows. 
It is easier to use the alternative expression for k as 

k{s,s') = p{s,s') - q*{s)r{s')q{s) . 
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Then, 



6k =6p-6{q*{s)r{s')q{s)) 

= p{s) X + - 6 {q*{s, s')p{s'Us, s')) 

= p{s) X + - 5q*{s)r{s')q{s) - q*{s)6r{s')q{s) - q* {s)r{s')6q{s) 
= p{s) X S(s) + 

- 5*{S) {p{s) - k{s, s')) - {p{s) - k{s, s'))5{s) - ^(s, s')^ {s')i* {s, s') 

= - f{s, s')^{s')^{s, s') - k{s, s') X S , (10.14) 



which is a direct analogy of (3.19). We have used the fact that for purely imaginary 
s we have, again 



S*{s)p{s) + p{s)5 = —5p + pS = p X S 



and 



s*(s)k(s, s') + k(s, s')s = — sk(s, s') + k(s, s')s = k(s, s') X S . 
Next, let us define the purely imaginary quaternion 

l*{s,s')5i{s,s') :=2Tis,s'). 

The real part of ^*{s, s')6^{s, s') vanishes since 3(s,s') is a unit quaternion. The last 
step is the variation with respect to T, computed as 

f{s,s')S^{s,s') = f{s,s')6[q*{s)q{s')] 

= f{s, s') [~q*{s)6q{s)i{s, s') + q*{s)q{s')q*{s')6q{s')] 



Thus, we find 



2T{s,s') ■.= lm{f{s,s')6i{s,s')) 

= -3*(s,s')S(s)3(s,s') + S(s') 
=: -Ad,.S(s) + S(s'). 



(10.15) 



Note the exact correspondence between this formula and (3.16) defining the variation 
Therefore, the variation with respect to 6k, gives 



^-^,Sk 

OK 




dU 
Ok 



s,s') X K{s,s')ds' , S(s))ds 



+ 




8U 8TI 

— (., s') - i{s, s)^{s', s)f{s, s')ds' , *(.) )d. . (10.16) 
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Analogously, 



//(.f,r.)d.d.^//(, 



3*(s,s')S(s)3(s,s') + S(s'))dsds' 
. dU 



Im 



dU 



ds',S(s))ds. (10.17) 



Collecting together the terms proportional to and ^{s) in the minimal action 
principle (10.7) gives the system (3.3), (3.4). The role of antisymmetric matrix Z{s, s') 
describing the nonlocal interactions in (3.4) is now played by the purely imaginary 
quaternion 



Z{s,s') = Im 



A Hamiltonian description closely following that of Section |5] can be developed in 
quaternionic form, as well. The cross products are then substituted by a correspond- 
ing product of the quaternions, with explicit formulas for the Lie-Poisson bracket 
closely resembling (5.3). Since the derivation is analogous to Section |5| it will be 
omitted from the exposition. 



11 Outlook for further studies 

This paper formulated the problem of strand dynamics for an arbitrary long-range 
intermolecular potential in the convective representation [TT] of exact geometric rod 
theory [12]. Its methods would also apply in the consideration of Lennard- Jones 
potentials and the constrained motion of non-self-interacting curves. 
After a quick derivation of the equations of motion by the Hamilton-Pontryagin ap- 
proach, the paper demonstrated and compared three different approaches to deriving 
the same continuum equations of motion for an elastic strand experiencing nonlo- 
cal (for example, electrostatic or Lennard- Jones) interactions. These were: (1) the 
Euler-Poincare approach; (2) the affine transformation approach and (3) the covari- 
ant Lagrange-Poincare formulation. The paper concentrated primarily on the case in 
which the strand is one-dimensional, which is the main object of interest for biological 
applications. However, these three approaches possess more significance and applica- 
bility than might be suggested by the one-dimensional developments illustrated here. 
For example, the geometrical considerations and nonlinear context of the present in- 
vestigation also apply in formulating the dynamics of the higher dimensional case. 
That is, when s has more than one component, the approaches discussed here still 
apply. 
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A change in dimensionality of s in equations (3.3) and (3.4) requires summing over all 
components of s-derivatives (instead of only the single s-derivative for the strand). 
Additional integrability conditions arise from the equality of cross-derivatives with 
respect to space and time that generalize equations (2.18) and (2.19). (In geometric 
terms, these are zero curvature conditions.) The extension to higher dimensions 
was discussed in the general setting treated in Section |2.4[ The higher dimensional 



options also figured in the covariant Lagrange-Poincare formulas (9.1), where div^. 
denotes derivative with respect to time and all dimensions of the space (taken to 
be one-dimensional in the paper). The extension to higher dimensions illuminates 
the geometry underlying the present one- dimensional case and may be expected to 
produce interesting applications in the dynamical description of biological membranes 
and other extended physical objects. While the equations take the same geometrical 
form in higher dimensions, their solutions will possess their own unique features. 

Besides passing to higher dimensions, future studies will consider both linear and 
nonlinear wave propagation on electrostatically charged strands, as well as the de- 
scription of nontrivial stationary states that arise from nonlocal interactions, such as 
for the VDF oligomers mentioned in the Introduction. 

Yet another interesting question for future studies concerns the possibility of enhanc- 
ing the internal structure of the rigid charge conformations. This will allow even richer 
dynamics than we considered here. While the resulting equations may be different 
(and more complex), the methods developed in this paper will still be applicable when 
the dynamics takes place in spaces that possess richer conformational structure than 
rigid rotations. 
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